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When you set out for Ithaka 
pray that your road's a long one, 
full of adventure, full of discovery. 
Laistrigonians, Cyclops, angry Poseidon 
don't be scared of them: 
you won't find things like that on your way 
as long as your thoughts are exalted, 
as long aS a rare excitement 
stirs your spirit and your body. 
Laistrigonians, Cyclops, wild Poseidon 
you won't encounter them 
unless you bring them along inside you, 
unless your soul raises them up in front of you. 


Pray that your road's a long one. 
May there be many a Summer morning when 
full of gratitude, full of -joy 
you come into harbors seen for the first time; 
may you stop at Phoenician trading centers . 
and buy fine things, 
mother of pearl and coral, amber and ebony, 
sensual perfumes of every kind, 
as many sensual perfumes as you can; 
May you visit numerous Egyptian cities 
to fill yourself with learning from the wise. 


Keep Ithaka always in mind. 

Arriving there is what you're destined for. 
But don't hurry the journey at all. 
Better if it goes on for years 
so you're old by the time you reach the island, 
wealthy with all you've gained on the way, 
not expecting Ithaka to make you rich. 
Ithaka gave you the marvelous journey. 
Without her you wouldn't have set out. 

She hasn't anything else to give. 


Andedfhyouefind*her poor), elthakacwonet thaveutooledniyou: 
Wise as you'll have become, and so experienced, 
you'll have understood by then what an Ithaka means. 


Kostas Kavafys 
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ABSTRACT 
The purpose of this study was to develop techniques for 
analysis of the hammer throwing and to analyse the throwing 
technique of world caliber hammer throwers. The development 
of the techniques included testing of different versions of 
the DLT 3-D cinematographic method and the derivation of 
formulas, based on rigid body dynamics, for analysis of the 
kinematic and kinetic parameters involved in rotation of the 
human body. 

The DLT was tested with and without mathematical models 
for image refinement and with calibration trees of different 
shape. A tree which was geometric in shape used with the 
basic DLT equations and with model IV for image refinement 
(Karara and Abdel Aziz, 1974) gave the best results of all 
the tests. The RMS error of the measured versus the 
Simulated coordinates of the control point were found to be 
0.24 -emenborithe Xsaxis; O0.d94emhfor ethecYraxis eand;40%26 cm 
for the Z-axis. The method was also tested for areas outside 
the calibration tree. It was found that a well constructed 
tree with image refinement models can be used for 
calibrating areas of larger volume than the one covered by 
the calibration tree. 

The hammer throwing data were collected during the 1982 
European Championship, with the best throw of the three 
medalists being analysed. Lagragian interpolation formulas 
were used for the time-match of the coordinates derived from 


the two films. The data were smoothed with digital filters. 
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It was found that the maximum velocity of the hammer 
occurred before the release point. None of the athletes 
achieved an optimal angle of release, while the closer to 
the optimal angle achieved by the winner could give him the 
gold medal with the same difference in throwing distance. 
The acceleration of the hammer was maximum in the last 
double support phase of the throw. It was the duration of 
the acceleration in this phase that was critical for a 
successful throw. Another parameter that characterised the 
analyzed athletes was their ability to increase the 
acceleration of the hammer during the single support phase 
of each turn or at least not to decrease it a great deal. 
This was accomplish by decreasing the moment of inertia of 
the body about the vertical axis. All athletes achieved 
faster single support phases than double support phases in 
all the turns. The latter was achieved by leaning backwards 
during the first half of the double support phase and by 
using the second half of this phase to initiate a fast 
rotation. The breaking of the horizontal movement in the 
last turn was critical for a successful throw. A relatively 
low center of mass of the body seemed to be appropriate for 
controlling the movement in all the turns. A tall and/or 
heavy athlete does not necessary have advantages against a 
shorter and lighter athlete. The length of the hammer should 


be shortened in order for athletes to perform mechanically 


better throws. 
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I. INTRODUCTION 


Biomechanics of sport is a discipline which studies the 

mechanics of human motion under the special condition of 

sports. Its content can be grouped into four major areas. 

1. Determination of the characteristics of an outstanding 
performer. 

2. Detection of differences between performers. 

3. Investigation of new techniques of performance. 

4. Validation and supply of data for mathematical models of 
performance. 

Biomechanists derive information about human motion 
mainly by means of: (a) Direct force 
measurements;(b) Electromyography; (c) Electrogoniometry; 
and (da) Cinematography. Direct attachments to the subject 
are required with methods (a), (b) and (c), while in (d) a 
Camera can be placed a considerable distance from the 
subject ensuring no physical interference. 

Cinematography has been a most useful tool for 
researchers in biomechanics and a most acceptable tool for 
subjects when data has to be collected in real life 
Situations. 

Most of human motion occurs in three-dimensional space 
rather than in a plane. However, film as a plane does not 
provide the analyst with information about the third 


dimension. Investigators in biomechanics have sought methods 
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of deriving three-dimensional coordinates of the body in 
motion by using films. In 1930, Bernstein introduced the use 
of a mirror in order to obtain a second image of the motion 
by using one camera. Since that time a variety of 
three-dimensional methods have been introduced. Within the 
last decade several investigators have drawn upon techniques 
from other disciplines in their attempts to resolve the 
problem. 

Photogrammetry is a discipline in which the fee of 
highly sophisticated metric cameras enable measurements 
based on the photogrammetric concept that the film, being a 
perfect plane, is a central projection of the object space. 
Assuming that the internal characteristics of the cameras 
are well known and by using a comparator and observations 
from different points, one can achieve spatial coordinates 
of a point in the object space. This is a normal 
photogrammetric procedure since the characteristics of the 
metric cameras are well known. 

In 1971, the photogrammetrists Abdel-Aziz and Karara 
presented a technique called "Direct Linear Transformation 
(DLT)" with which one can achieve photogrammetric 
measurements by using non-metric cameras. This is a two-pass 
regression technique; a set of control points (calibration 
tree) in the object space is used to calculate the external 
and internal characteristics of the cameras. The calculated 
coefficients are used together with digitized coordinates of 


the points of interest to achieve the spatial coordinates of 
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the points. 

Shapiro (1978), Walton (1981) and Miller et al (1980) 
tested the technique for biomechanical research and found 
the resulting measurements to be highly satisfactory. The 
calibration of the cameras in these studies was achieved by 
using a three dimensional calibration tree. In 1980, 
Woltring presented a technique in which calibration could be 
achieved by filming a plane grid with control points. The 
grid was rotated to different positions and filmed. These 
different positions of the grid offered a 3-D calibration 
tree. This method has the advantage of ease of construction, 
but it is difficult to have a large grid in rotation without 
deformation. 

The relationship between the volume of the calibration 
tree and of the calibrated area has been of main concern 
among the researchers involved with DLT. To resolve this 
problem Dapena et al, (1981) presented a version of the DLT 
with which one can use a control object of unknown shape to 
calibrate large areas. Weak points of this technique 
included complexity in the calibration process, and failure 
to account for refinement of the image caused by lens 
distortion and film deformations. It is generally accepted 
that data collected with non-metric cameras encompasses 
errors caused by these factors (Shapiro 1978, Walton 1981 
and Karara and Abdel-Aziz 1974). In 1974, Karara and 
Abdel-Aziz presented six mathematical models which when used 


together with the basic DLT model account for image 
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refinement. These models have not been utilized in previous 
biomechanics research. 

Assuming that the body is a mechanical link system and 
if the coordinates of points of interest of the body can be 
accurately determined,in. 3sD space, one.can.use classical 
mechanics methods to investigate human motion. 

Hammer throwing originally appeared as a hunting 
method, later as means of waging war asS a weapon and 
finally, in todays "civilized world", as an Olympic event. 

Despite its change in purpose, the movement pattern has 
remained basically the same. A heavy object mounted at the 
end of a stick, or a belt, or a string, is released after a 
few rotations in order to kill a dinosaur, to wound a 
Goliath, or to break the world record. According to Howard 
Payne, 

"Hammer throwing as we know it today, has had a 
short history if we date it from 1687 whemethe 
Americans drew up rules setting the circle 
diameter at 7 feet., the hammer length at 4 ft. 
and the weight of the ball, chain and handle of 
the implement at 16 1b. However, if we include 
sledge hammers, shafted implements and wheel 
hubs, then we can go way back to trace the 


Origins -some historians think to 2000 B.C. when 
at the Tailteann Games in Ireland the "Roth 


Cleas" or "wheel feat" was contested, ...", 
(Payne 1969:9) 


The specifications of today's Olympic event, according 
to Rules and By-Laws of the Canadian Track and Field 
Association, are presented below. 

The hammer is composed of a sphere (the head) and a 
grip connected by a steel wire, as in Figure 1. The minimum 


weight of the hammer is 7.257 kg. and the complete length 
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must be between 117.5 and 121.5 cm. A glove can be worn by 
the athlete to protect his fingers while throwing. 

The Head, must be spherical with a diameter between 
102mm. and 130mm. It should be made from solid iron, or 
brass, Or any metal not softer than brass, or a shell of 
such metal filled with lead or other solid materials. If a 
filling is used it must be inserted in such a manner that it 
is immovable and the distance of the centre of gravity from 
the centre of the sphere is less than 6mm. 

The Grip must be triangular and may be either of single 
or double loop construction, but must be rigid and without 
hinging of any kind and so made that it cannot stretch 


appreciably while being thrown. 


Figure 1: Dimensions of the Hammer 


The Connection Wire must be a single unbroken and 


Straight length of spring steel wire with a diameter no less 
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than 3mm. It may be looped at one or both ends as a means of 
attachment. The connections of the wire with the grip and 
the head shall be such that they do not increase the length 
of the hammer while it is being thrown. 

The throwing circle must be made from concrete and have 
a slightly rough surface. The diameter of the circle must be 
2.134 meters, andthe throwing arc 40°. For protection 
purposes, a cage with certain specifications is required. 
Figure 2 1s a schematic representation of the throwing 
circle and the protection cage. 

The following is a description of a right-handed 
throwing technique. The athlete stands at the back of the 
Circle with his feet almost parallel and facing the opposite 
direction of the throw. The gloved left hand is placed in 
the handle so that the grip rests along the middle phalanges 
of the fingers. The right hand is placed under the left. The 
athlete swings the hammer a few times around his head while 
his feet remain almost stationary. At the end of the last 
swing the thrower initiates a whole body turning action with 
hic feet, rotatapg with his weight over hiswbett foor. The 
gradually accelerated rotation carries him to the front of 
the circle after three or four turns and he uses the last 
turn to give the hammer a final large release velocity. The 
plane of the orbit of the hammer throughout is inclined at 


an angle so that an optimal angle of release can be 


achieved. 
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Figure 2: Hammer Throwing Circle 
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The hammer is a projectile and assuming negligible 
aerodynamic factors or Coriolis forces duringeidtoht,, the 
final distance of a throw is given by the following 


equation. 
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where: L = Horizontal distance of the throw. 
U = Velocity of release. 
6 = Angle of release. 
h = Height of release. 
g = Acceleration due to gravity. 


Also, by differentiating (1) with respect to 6, with 
constant velocity and height of release the distance L will 


be maximum if: 


yy? 

sin?@ = —————————_- (2) 
GOS eh os) 

Thus, assuming an optimal angle of release it is reasonable 

to consider the velocity of release as the most important 

factor affecting the final distance of a throw. 

The world record for hammer throwing has shown greater 
improvement than the records for any other field event. This 
is mainly a result of continued improvement in technique. 


The world records in this event and in the discus throwing 


event for the last four decades are shown im Fabletd. 
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TABLE 1 


WORLD RECORD OF HAMMER AND DISCUS THROWING 
POR 1950370 1980 


Year Hammer Discus 
1950 59.80m. 57 <o3m. 
1960 J0.33m 60.46m. 
1970 75.84m. 68.40m. 
1980 81.80m. UM ems 


Although the technique has been continuously evolving, 


there is limited literature available concerning the event. 


Publications which are available are almost entirely 


Qualitative in nature. The lack of quantitative analysis may 


be a result of one or both of the following. 


pe 


The lack of interest in the event by non-Europeans. In 
the last two decades, the event has been dominated by 
athletes from the U.S.S.R. The three medals in this event 
in the last three Olympic Games have been collected by 
athletes of that country. 

Hammer throwing, a complex event, is basically a 
rotational movement, therefore, the value of a 
two-dimensional film analysis is questionable. The 
physical space required for this event increases its 
complexity. It is obvious that a controlled laboratory 
environment is not the best place for throwing the 
hammer. 


The underlying responsibility of the sports 


biomechanist is to provide the coach and the athlete with 


c re) “Go pm? 
o VSS 


evoeit 
mee. $2 i008 BE 


8%, 06 mee, OF vs 
wob8s eae 
mstcts ery tan: 2 

é A 


png iviove, yfaueumt tage mged 2cH aiptarags ons dl 

.tnsvs9 edo eninisonog ay iantine ‘pagiateeatt a “ 
yleur ins, deoiLs S75 2idatiata ade apie 

ysm efeylens oytiedioasup Ho: Hoaqt adP ee 

} eridvotito ss ott “edsbaine/ 


S9qo1vS-nen ed jneve sty bi: 


ed 
“a 


yd betedimos yaaa See Broval ena: 
| ps i 
tnev> erde ne 2 eben soutabadts phe 


yd Hasosllos feed SK ath a 9 2aigi Oe 


yo ‘se 


& to ast ca ‘Bdz re { 4 peat 
ai? .sldenotieayp at ie . : obs 


ati es25e tant inate eats, ’ 
g¢roder0dsk bsiiornod’s salah a 38 seats 
ed¥ paiwoxdd 303 sosiqteud ez 20% ion ab trsmeesivn 


r Siaeqa siz to yi 
fitiw etelite- edd bus | 


10 


useful information about kinematic and kinetic parameters of 
the skill. The lack of such information about hammer 


throwing performance inspired this study. 


A. Purpose of the study 

The purpose of this study was to develop a method for 
analysis of hammer throwing and to analyze the technique of 
world caliber hammer throwers as it is performed in a 
competition. In order to develop the method for analysis, 
different versions of the DLT 3-D cinematographic method 
were tested. These tests included the basic DLT method with 
different mathematical models for image refinement, as well 
as different calibration trees for the derivation of the 
calibration coefficients. Rigid body dynamics were used for 
the derivation of the formulas to investigate the parameters 
for analysis. Temporal, kinematic, and kinetic parameters of 
(a) the hammer, (b) the athlete, and (c) the system, athlete 
plus hammer, were analysed. 

A comparison between between the different subjects was 


made based on the analysed parameters. 


B. Limitations 
The conduct of the study was limited by the ability of 
the researcher to determine the location of selected body 


and object points, and also by the number of subjects. 
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C. Delimitations 

The study was delimited to the three medalists of the 
1982 European Championship, to a filming frame rate of 100 
frames per second, to the motion of the hammer and the 


System after the preliminary swings. 


D. Definition of terms 


Calibration of an image is the process of determining a 
number of unique coefficients which when substituted 
into an equation or a set of equations define the 
object to image transformation. 


Calibration coefficients are the unique numbers which when 
substituted into equations define the object to image 
transformation. 


Calibration tree is a structure which provides a set of 
points (control points) with well known spatial 
coordinates and they are used for the calculation of 
the calibration coefficients. 


Central projection is a mapping of points in which each 
point, its projection and the projection center are 


collinear. 


Comparator is a precise optical instrument with viewing 
optics for obtaining coordinate pairs from a 


secondary image. 


Digital filter is a set of algorithms which form an 
analytical tool for removing or isolating certain 
frequency components from a digital signal. 


Digitizer coordinates are a pair of coordinates which 
describe the location of a point with respect to 


digitizer reference frame. 


Digitizing board is a precise instrument which is used to 
obtain coordinate pairs of points froma film. 


Direction cosines are a set of numbers which define the 
orientation of a vector in an orthogonal system of 


axes. 


Direct Linear Transformation (DLT) is a two-pass linear 
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regression technique with which three dimensional 
coordinates of a point in space can be achieved by 
using films from different observations. 


Double support is the phase of the throw when the athlete 
has both feet on ground. 


Film deformation is the image deformation caused by changes 
in the geometry of the film. For example the lack of 
a perfect plane when the film is in the camera or 
DPFrOR eC COT 


Focal center is the point of the lens through which all the 
optical rays appear to pass. 


Instantaneous Center of curvature is the center of the 
instantaneous circle described by the orbit of a 
point which traces a rotational motion. 


Interpolation is a numerical method in which a function is 
evaluated at some point within a specified interval 
of the independent variable based on neighboring 
values of the dependent variable. 


Least square technique is a numerical method in which the 
parameters of a mathematical model are determined by 
minimizing the sum of the square differences between 
the observed and modeled values of the dependent 
variable. 


Lens distortion is the geometric distortion of an image 
caused by the imperfection of the lens. 


Metric camera is a precision camera for recording 
photographic images in photogrammetry. Some of the 
characteristics of this camera are: A fixed focus; 
Symmetric lens with minimal optical distortion; The 
ability to hold the film flat; A diaphragm shutter 
located within the lens; A precisely calibrated inner 
chamber for which the parameters of interior 
orientation are known. 


Newton’s iterative method is a numerical analysis method to 
find the roots of an overdetermined system of 


non-linear equations. 


Newton’s method is a numerical analysis method for 
determining the roots of an equation. 


Non-metric camera is a camera which does not meet the 
requirements of a metric camera.(see metric camera). 


Overdetermined system of equations is a system of equations 
in which the number of equations exceeds the number 
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of unknowns, 


Photogrammetry is science of obtaining measurements from 
frims. 


Precessional motion is the motion of a spinning top where, 
as the top rotates about its axis of rotation, the 
axis of rotation rotates about another axis 
(vertical). Nutation is the special case of 
precession when the angle of precession does not 
remain constant. 


Radii of curvature is the distance from the point which 
traces a rotational motion to the instantaneous 
center of curvature. 


Random error is an error which is created by the 
inconsistencies of the measurement process. 


Single support is the phase of the turn in which the athlete 
has only one foot on the ground. 


Spatial coordinates are a pair or triplet of numbers which 
define the location of a point with respect to a 
reference frame. 


Systematic error is an error which is consistent in the data 
and caused from the malfunction of an instrument or a 
mathematical model. 
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II. REVIEW OF RELATED LITERATURE. 


The review of literature related to this work is presented 
in the following sequence: 
2.1 Three-dimensional cinematographic techniques. 
2.2 Smoothing techniques. 


2.3 Hammer throwing. 


A. Three-Dimensional Cinematographic Techniques. 
The techniques found in the literature can be divided 
into four basic groups: 
1. Single camera technique. 
2. Mirror technique. 
3. Stereometric technique. 


4, Multi-camera techniques. 


Single Camera Technique. 

Plagenhoef (1968) presented a method that utilized the 
real length of body-segments, the camera-to-plane of motion 
distance and the length of the segments as they appeared in 
the film in order to calculate the spatial coordinates of 
points by using trigonometric functions. The same technique 
has been used by Riley et al (1978). In their study, a 
simulated computer graphics method was used to determine the 
third coordinate. The techniques described above give a 


rough estimation of the third coordinate. In analysis of 
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Single segment movements, the results can be reasonable, but 
the lack of accuracy is obvious in more complicated 


analyses. 


Mirror Technique. 

This sophisticated technique was introduced by the 
Russian physiologist Bernstein (1930). A plane mirror is 
placed in such a way that the object space plus its image in 
the mirror are recorded on the film. With this technique, a 
researcher can obtain two images from different points of 
view by usSing only one camera. Although the technique has 
the advantage of absolute synchronization of the two 
"cameras", it suffers from the disadvantages of: (a) small 
image, (b) mirror surface distortion and (c) it can be used 


only in a laboratory environment. 


Stereometric Technique. 

This technique is based on the same principles as those 
of binocular vision. The eyes, being approximatively 7 cm 
apart record an object from different angles. The two 
different images are reconstructed in the brain giving the 
impression of the space more than of a plane. In practice, 
two photogrammetric cameras are placed side by side with 
their optical axes parallel to one another. For the 
calculation of the X,Y,Z coordinates, standard equations 
have been derived. Ayoud et al (1970) were the first to 


Suggest the application of this technique for obtaining 
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three dimensional spatial coordinates of human motion. 

This technique requires photogrammetric cameras with 
well known internal characteristics. Such cameras are very 
expensive, and consequently the technique appears to be 
unrealistic for the standards of a biomechanics laboratory. 
Moreover, the fact that they are still cameras and that an 
external device is necessary in order to obtain a multiple 
image (e.g. stroboscope) limits the technique to laboratory 


environments. 


Multiple Cameras Technique. 

This technique employs more than one camera to record 
the motion, and different algorithms can be found in the 
literature to reconstruct the three dimensional coordinates 
of the body in motion. 

Noble and Kelley (1969) used three cameras to determine 
the three dimensional coordinates of a moving ball 
describing the path of a right circular helix. Two cameras 
were placed in the horizontal plane, 90° out of phase with 
one another. The third camera was positioned directly above 
the apparatus. The timing of the films was synchronized by 
firing a flash bulb at the beginning of the filming. Each 
film gave two coordinates and by using different scale 
factors, the X,Y,Z coordinates were reconstructed. However, 
the errors presented by the authors (35% -0f the criterion 
value in acceleration) deem the appropriateness of this 


technique within biomechanical research questionable. 
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Anderson (1970) presented a method which utilized three 
Cameras, one sighted along each of the three orthogonal axes 
X,Y,Z. The two cameras in the horizontal plane were used to 
calculate the X,Y,Z coordinates while the camera in the 
vertical axis was used to provide correction factors for 
perspective errors. This method necessitates the recording 
of each point by all three cameras. 

Miller (1971) presented another technique which also 
utilized three cameras. The cameras were not placed on the 
three orthogonal axes, but instead, in a horizontal plane 
and 120° out of phase with each other. The intersection of 
the optical axes of the three cameras was used as the 
Origin. Three targets were carefully positioned beyond the 
Origin and aligned with the optical axis of each camera, 
providing reference coordinates of the origin. The space 
coordinates of the cameras and the distance between the lens 
nodal plane and the film plane inside the camera are 
required to use the algorithm for the reconstruction of the 
X,Y,Z coordinates of the filmed object. This technique has a 
general applicability and provides an exact analytical 
solution to the perspective problem. However, because it 
requires a specific placement of the cameras and of the 
targets, it is difficult to use in real life events and 
especially in a competition situation. 

Bergemann (1974) emphasized the importance of camera 
placement to the production of accurate spatial coordinates. 


The two cameras were placed in the horizontal plane with 


ae As 
of E = 
fi2 


eetds basi lite ott podjom 5 , tasndate 

eoxs Ensmeetete agsdd ade 19 do6e 1 

ost been or3W snes lg fegngsi vod -oft © i 

ada nt BISMBD. eft, elite sssentelgele sai 

102 Budd 962" noi?5e7707 Shi vorg o3 bpau 2502 

on ih103e7 sat es¢siteasasr Gotta ene proeie 

sh es: eres, asi die oe 

ate igndw eunsriine Isiigg8" ae en 

git 0 ssaete Jon SEX asisNisa, ant. apr nue 

atela ([a4tdosized ant bate 2 tee ical (Spe a 
rere 1Ss9n2 (snk Peet gaia en asst 
as 


éf3 an Beeus ebw, estsins> 965d" eid to eens =] 


i Secu 


aig ono ay ad pemokttecs oP pie ynin 
, 679R59 seo gied 20 wine igatigot rey é, 


st 


soeda sit aie sai to session 


sit lo: norssutte OIA Pies ba ae rou 
s.26d sup indoes hi? aoetdo! ‘pamits 

isslaylens: tDaxe 16 scsi seh 

4% sausosd ‘ee, ; oobsory mesione se oti 

6f4 to brs seuhiss fats So siemens iq’ = 

bis: efaeve stil [est ni sit ba eguats9n6 ei Ye st 

nos aude ined dé seqnos at yliekesqes 

Gtesso lo-soned20Qmi oft nesiesiiqns «4 (A721) onsmegist eo 

waessriwo? larssqe stenvose to nel s2dbo3a adz od Ce) ce | 


: 
ase 7 = 
2 


¢ Seidl 
+0 
- ion ~ 


re 


18 


their optical axes intersecting at a common Origin. For the 
placement of the cameras, standard surveying equipment was 
used. Equations were derived to calculate the position of 
several arbitrary points on a coordinate grid. A similar 
method was presented by Van Gheluwe (1974). Pensore et al 
(1976) have shown that it is no longer necessary for the 
optical axes to intersect, although a theodolite is required 
to obtain precise spatial information for the cameras and 
the reference points. 

Van Gheluwe (1978) reported a technique involving 
camera placement in any position relative to the subject 
being filmed. The basic principle of the technique relies 
upon the explicit mathematical reconstruction of the 
position of the cameras in the space using the known 
life-sized coordinates and the image-coordinates of certain 
reference points. A steel tri-axial reference frame was 
placed in the filmed space to provide appropriate 
information for the calculation of the position and for 
orientation of the cameras. The accuracy of the presented 
results were within the acceptable limits. However, it is 
not known whether the points chosen for the reliability test 
were other than control points. Morever the method does not 
account for image refinement from lens distortion or film 
deformation errors. 

In 1971, Abdel-Aziz and Karara presented a method that 
allowed the use of non-metric cameras (cameras with unknown 


internal characteristics) in stereo-photogrammetry. This 
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method is the "Direct Linear Transformation" (DLT) method. 
The position and orientation of the cameras were derived 
through the use of control points in the space instead of a 
comparator. A minimum of six non-coplanar control points was 
required for the solution. Six mathematical models were 
presented by the same investigators (Karara et al, 1974) to 
refine the image from lens distortion and film deformation 
in non-metric photographs. 

Shapiro (1978) tested the DLT method with high-speed 
cinematography. A "tree" with 48 control points was used to 
calibrate the cameras. It was found that the unknown spatial 
coordinates of a point can be located with an average error 
Sf1t0;4emafor+the Xsaxisp +0 s5emefor-the-Y-axis+and.+0.5em 
for the Z-axis. In a dynamic test the acceleration of a ball 
inetree-falluwas found: to-yieldsybetween-79.5,tops10em/s’. 

In an extended investigation of the three dimensional 
problem, Walton (1981) arrived at the conclusion that the 
precision with which the object points were located with the 
DLT method, was adequate for most analyses of human motion. 
He suggested that accuracy of the measurements can be 
increased by using a redundant number of control points and, 
if possible, more than two cameras. 

Miller et al (1980) presented the use of the DLT method 
together with a method for obtaining spatial kinematic 
Parameters of segments of biomechanical systems. 

The problem of producing accurate three dimensional 


coordinates has been of continuing interest. Virtually all 
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of the above studies validated their procedures by filming 
an artificial frame. Walton (1981), in one of his validation 
procedures, used the body landmarks of a trampolinist to 
calculate the acceleration of the C of M of the athlete due 
to gravity. He found that the largest error in any of the 
computed values for "g" was 1.58 m/s?. He also found that 
the mathematical model of the DLT solution with the image 
refinement model increased the accuracy of the measurements. 
He suggested that the data obtained through cinematography 
Should be filtered to reduce inherent errors associated with 
film analysis. 

Woltring (1980) tested the same technique. For a 
calibration tree, he used a 2-D plane with control points. 
The plane was rotated about an axis and was filmed in some 
specific angular positions. This technique can be very 
useful ina laboratory enviroment or for calibrating small 
areas ypebutvit’ is®dif£fteult*to* find = construetion” material 
which will not undergo deformation when the plane is of a 
large size and must be rotated. 

A alteration of the DLT method was presented by Dapena 
et al (1981). With this method the use of a calibration tree 
is no longer necessary and the calibration of the cameras is 
achieved by the use of two calibration crosses and a set of 
vertical poles which are placed in the filming area. The 
authors validated this method by using a 24-point grid. The 
average Root Mean Square error of the measured versus the 


Simulated coordinates of 15 points of the grid was 1.50cm 
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for the X-axis, 1.30cm for the Y-axis and 0.60cm for the 
Z-axis. This model did not account for the image refinement 
components and in a large filmed area the image deformations 
from the lens or the film itself would be expected to be 
large. Further, these investigators did not justify the 
deletion of 9 points of the tested grid (15 points of the 


grid out of 24 were used for the validation). 


B. Smoothing Techniques. 

Since there is always a certain amount of error 
associated with film analysis, it is necessary to use 
numerical methods which reduce this kind of error. There are 
several different methods of data smoothing utilized in 
biomechanics research. Early researchers used manual 
smoothing methods to draw smooth curves through successive 
data points (Miller and Nelson, 1973). Finite differences 
have been used to calculate first and second time 
derivatives of displacement data (Miller and Nelson, 1973). 
Plagenhoef (1973) used a polynomial approximation method to 
Smooth film data. Using displacement data, he determined the 
coefficients of a polynomial which best fit the data points. 
The order of the polynomial could be altered depending on 
the regularity of displacement data. The polynomial was 
differentiated to obtain estimates of velocity and 
acceleration data. 

Reinsch (1967) presented the use of spline functions to 


smooth experimental data. The use of cubic splines has been 
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widely accepted by biomechanics researchers and it has been 
reported that the acceleration values calculated by this 
method for film displacement data were more accurate than 
those derived from other methods in the past (Zernicke et 
al, 1975; MacLauglin et al, 1976). To solve the problem of 
acceleration of the displacement-time data becoming zero at 
the knot points, Zernicke et al, (1976) suggested that three 
extra points at the beginning and at the end of each data 
set should be used. 

The estimation of the degree of smoothing is an 
additional problem for the researcher. MacLaughlin et al 
(1976) suggested that the calculated average error in 
measurement of distances from film can be used for the 
estimation of the degree of smoothing. Further, Zernicke et 
al (1976) suggested the use of the resultant acceleration 
curve to alter the degree of smoothing. 

A more recent smoothing method was reported by Pezzack 
et al (1977) and it is known as the digital filter method. 
They used a second order Butterworth low pass recursive 
filter to reduce the noise inherent in the signal of film 
data. The acceleration curve of a freely rotating arm was 
used to validate the accuracy of smoothing with this method. 
Finite differences and Chebyshev least squares polynomials 
were tested in the same experiment. An accelerometer was 
connected to the arm to provide a criterion acceleration 
curve. It was found that the digital filter followed by 


finite difference differentiation was the only method to 
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accurately reproduce the acceleration time curve recorded by 
the accelerometer. Spline functions were not tested in the 
above study. However in their discussion, the authors 
doubted that spline functions would give more accurate 
results than Chebyshev polynomials. 

The aforementioned studies present methods for 
smoothing data in biomechanical research. However the 
question of whether a signal of certain frequency is a data 
Signal or a noise signal remains unanswered and the 
estimation of the error inherent in data analysis remains a 


subjective consideration. 


C. Hammer Throwing. 

The first attempt to biomechanically analyse hammer 
throwing technique was by Lapp (1935). He found that the 
hammer accelerates when the athlete is in the double support 
phase and decelerates during the single support phase. 

Samozvetov (1974), in a study of the acceleration of 
the hammer found that the duration of acceleration of the 
hammer was critical for a maximum distance throw. He also 
found that the acceleration phase corresponds to the double 
Support phase of the thrower. 

Kuznyetsov (1974) used wire tensors attached to the 
handle of the hammer to measure centripetal forces caused by 
the effort of the thrower to propel the hammer. This was 


used as a strength measurement method but results were not 


presented. 
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Dyson, (1977) claimed that, throwers can accelerate the 
hammer in the single support phase by pulling it downward. 
"Vertical downward acceleration is achieved by permitting 
the body weight to drop just before the hammer head itself 
drops in its swings or turns. Here, therefore, the limiting 
factor to acceleration is the thrower’s weight - which can 
be lowered with or without contact with the ground." (Dyson, 
1977). If this statement is correct then the weight of the 
athlete is another critical factor in hammer throwing. 

Elmendorf (1978) studied the relationship of turn time 
and distance thrown. He found that turn time and distance do 
not have significant relationships in a three turn throw, 
but they are significantly related when the athlete uses 
Pour turns. 

Ariel (1978) compared kinematic parameters of the six 
best throwers of the 1976 Olympic Games with those of an 
average group. He found that the average group had a longer 
turn time for each turn. He also found that the double 
Support time of the Olympians was longer that those of the 
second group. 

Dapena (1981) suggested that the tangential component 
of the wire pull force, rather than the tangential component 
of weight of the hammer force is mainly responsible for 
changes in the velocity of the hammer. He found that for a 


throw of 67.50 metres, the maximum total wire pull was 2750 


Newtons. 
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The above studies have presented limited quantitative 
information about the hammer throwing technique. The need 


for further study is evident. 


Ill. METHODS AND PROCEDURES 


The methods and procedures utilized in the present Study are 
described below. The DLT method was tested together with 
different calibration trees and the six mathematical models 
for image refinement presented by Karara and Abdel-Aziz in 
1974. The model which gave the best results was used for the 
data collection and data reduction of hammer throwing. The 
digital filter method was used for the reduction of the 
noise (Pezzack et al (1977) and Winter (1979)). Rigid body 
dynamics methods were employed for the analysis of the 


hammer throwing technique. 


A. DLT Method 

The DLT is based on the basic theoretical concept of 
photogrammetry: the photograph, being a perfect plane, is a 
central projection of the object space. Analytical 
representation of the geometry of the technique can be found 
in Appendix A. It is normal photogrammetric practice to 
determine the transformation of coordinates from the 
secondary image to real life coordinates by careful 
calibration of the metric cameras and comparators. However, 
when non-metric cameras are used, an external calibration is 
required. For this purpose, a set of known points named 
control points (CP) are located in the object-space. If X, 


Y, Z are the spatial coordinates of a point in space and U, 
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V its digitized coordinates, the transformation is obtained 
by utilizing the following equations named basic DLT 


equations: 


L,X + b3¥ + L32Z + Ly 


(3) 
Bex thebmg Vino) alate’ 


Get PobeYeteL7Z #4 b3 


Vv + AV = SS (4) 
Lexeseloatrt LBanSbet a 


where; X,Y,Z are the coordinates of the point in space. 
U,V are the digitized coordinates of the point. 
L,(i=1,...,11) are coefficients named transformation 
coefficients. 
AU,AV are image refinement components of lens 
distortion and film deformation. 


A minimum number of six CP is required for the solution 
of a system of twelve equations and eleven unknowns. For 
better accuracy a redundant number of control points is 
advised. In this case an over-determined system of equations 
has to be solved and a least Square method can be used. In 
order to calculate 3-D coordinates of a point, at least two 
observations (cameras) are needed. 

Once the transformation coefficients of the two cameras 
have been determined and the digitized coordinates of a 
point in the object-space from both cameras are given, the 
DLT equations are used for the determination of the spatial 
coordinates of the point. 

Six mathematical models have been given by Karara and 


Abdel-Aziz (1974) for the image refinement components 
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(AU,AV). These models can be found in Appendix A. The first 
model does not change the form of the DLT equations, because 
the eleven basic coefficients account for the linear 
components of lens distortion and film deformation. When one 
of the left five models is used, the number of calibration 
coefficients and the minimum number of CP are increased as 
shown in Table 1, Appendix A. With utilization the models II 
to VI and with a redundant number of control points, the 
system of equations to be solved becomes an overdetermined 


system of non-linear equations. 


Validation of DLT Method 

The DLT method was tested by utilizing three different 
calibration trees and the five mathematical models for the 
image refinement. Two Photo-Sonic IPL phase-locked cameras 
were used for the filming in all the tests. An Angenieux 
12-120 zoom lens was mounted on each camera. 

A computer program was written in HPL language to 
calculate the DLT calibration coefficients and to simulate 
the spatial coordinates of the CP. The DLT basic equations 
without any image refinement component were used as a first 
step. The least square method presented by Scheid (1968, 
pp.375) was used to solve the overdetermined system of 
equations. As a second step, the five mathematical models 
for the image refinement were sequentially introduced to the 
DLT equations. The digitized coordinates of the center of 


the projected images were calculated and were used as the 
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point of symmetry of the frames. Based on these coordinates 
and the digitizer coordinates of the CP the position vector 
of each point was defined with respect to the point of 
Symmetry. The overdetermined system of non-linear equations 
was solved by expressing the equations in truncated Taylor 
series and then by using Newton's iterative least squares 
method. The transformation coefficients derived in the first 
step were entered as initial values for the unknowns at the 
first iteration. In the first iteration the image refinement 
coefficients were assigned a value of zero. The computer 
program was set to attenuate the iterations when the sum of 
the absolute values of the residuals was smaller or equal to 
Bhetctiteria vaduerof)] 104182 e4ZRP Sod 0ehe) soorowhenathe 
system stopped converging. With the models II, III and IV 
the system satisfied the first criteria. With model V the 
system stopped converging when the criteria value was 

=R < 10-° and with model VI the system stopped converging 
when the value was IR < 10°%. This might be a result of the 
least square method used for the solution and of the 
precision of the computer (over 18 unknowns, 68 equations 
and requested high precision). 

The first tested tree was a 2-D frame with dimensions 
2X2.5 m (Woltring, 1980). The tree was constructed of metal 
tubes. table-tennis balls were used as control points and 
the tree was filmed in different angular positions. The 
Cameras were placed approximetaly 20m from the frame and 10m 


from each other. This particular disposition of the cameras 


es 


ah 
- fh 


astanthbroe> sega, nO beget a 


teioev neitiaog of? Wvsds to ge at 
20 dnfod: ‘Sd2 0 tages soe 


; ie ; 
enotwteups tésniic mort 10 motege, ‘pee ve 


yolvseT pstesaust. hi eno Fiaups. eae ( ia & fy 
pin 
asusupe teases avideveds & Gene 


, (yi > 
avid ed3 af bavi: “Sue tor hisos aod asqse 19, shit, 
add 3s. zavondal SAF 302, eer by ‘tap siqi aa & 

oi Paan 


inamanits: seem aha norsHisse ‘aks aa 


tsiuqnad enT 0795 to° St fey 6) oesente ar 

26 mue ond nedw  enorse gisti sg% » stiytadanie4 Jen 2s iw 
oF Laups ito roltsma. 2esawv bsuicenaaeh ae Torani a ub 
2) ete, to aula i 

4) >a 


c , 7s P . ‘ aa) 
" ui v rr A . ’ “ a5, yy . ‘ 4 “I 
vi. big D01h, 21 eepGom og? tate vpdepdauitob 


® 


‘ets astwo a .,¢*'*0! 


oe 


Ney 


ere A ;| aot ee) 
andy VY Debomatin® .sttsstiovre ie aac 
: \ 


sista botWota matte eal ie | 
sit to. tides) = Sd digta Baa 7 be AS eule 
end 16. sic eaten Per re 4 ge | 
ghoktveups Bd. ~eqionian 8f raved ‘aaduqines ~ 

| Aagiaiesera: be 
ghotenemib apie emani a-% 5 sewi9073 Beda 
Lesem. 96 Sedoursenos es* send: anit /CORET mia 
bas-atnieg Lottnes cs Seay. eraw atied etfnes-ediiad 
sit a 36 lupas thibanig oer nb) bamtia gow Band a 

set pun aaekostaimgll 05: ~Leskeitnnoi eit temalaae 
gaxemnso’ a¢3 to 7 be sasuottted eidt.% narnia 


30 


is not a critical factor in the accuracy of the results 
Walton; 1981; Miller let aly 1980) ¢itThe scamerascwéreosethto 
operate at a speed of 100 frames per second (fps). The 
measurement of the mean error was approximately 10 cm. This 
large error was probably caused by the deformation of the 
tree.sThese results indicated*that, for calibration of large 
areas and with non-expensive materials, the use Of Tan2=D 
tree wasS not viable. 

In March 1982 a new calibration tree, a steel frame 
with dimensions 2.5 x 2 x 2 m was tested. Thirty-five 
randomly distributed table-tennis balls were mounted within 
the frame by means of strings. The reference frame and a 
golf ball in free fall released from the top of the frame 
were filmed. 

The overall Root Mean Square (RMS) differences of the 
X,Y,Z coordinates simulated by the DLT method with the five 
mathematical models and the measured coordinates were 
calculated. Refinement with model III gave slightly better 
results than the other models. (RMS(X) = 1.12 cm., 

RMS (¥)©=01.09 cm:, RMS(Z) = 0.81 cm.) 5 All othe *models gave 
more accurate results for the Z-axis (vertical) than for the 
Other axes. This might be a result of the perspective errors 
of the two cameras which affected the X and Y axes or the 
precision of the instruments used for measuring the 
coordinates of the points in the experimental session. 

A golf ball in free fall was digitized in every frame 


of both films and the DLT equations together with Model ITI 
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for image refinement were used to simulate the Xa oematicerz 
coordinates of the ball. The time-displacement data of the 
ball were filtered with a second order double pass 
Butterworth recursive filter with a cut-off frequency of 8 
Hz. A first central difference method, (Miller et al 1973) 
was used to calculate the acceleration of the ball. This was 
Found to be between -9.5 m/s? and -10.2 m/s?. (Criterion 
Matte Of Quavity = -9.81 m/s?). 

The results of the above test were complementary to 
results of Similar studies and were, therefore, 
Satisfactory. Several decisions were made based on the above 
test. To improve the accuracy of the results derived with 
the DLT method, the cameras must be placed as far as 
possible from the object space to avoid perspective error. A 
pre-built calibration "tree" with precisely measured control 
points must be used. For better precision and faster 
measuring, the calibration tree should be geometric in 
shape. This can provide some CP which do not require actual 
measurement since the spatial coordinates can be calculated 
with respect to the coordinates of other points. 

The third tested tree was the calibration tree used for 
the data collection of the hammer throwing and was an 
aluminium frame of 3x3x3 m (Figure 3). It was made from 12 
aluminium sticks with a square cross-section (2.5x2.5 cm.). 
The sticks were parallel or perpendicular to the horizontal. 
The aluminium sticks were chosen instead of other material 


because they are easy to transport, easy to work with, and 
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Figure 3: Data Collection Calibration Tree 
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because a three meter stick, when mounted by two points one 
meter from each end, presents no deformation. Rigid L-shape 
corners were used to improve the stability of the 
connections. The connections and the end points of the 
Sticks provided a set of 33 CP plus one which was used as 
the space fixed origin. The coordinates of these points are 
presented in Appendix D. First, all the points were used to 
calculate the coefficients with each model and the X, Y, Z 
coordinates of all the points and the location of both 
cameras relative to the space fixed origin were simulated. 
Second, 27 points were treated as control points and 6 
points (points 28 to 33 in Figure 3) were treated as 
unknowns. The calibration coefficients were recalculated 
with all the models and the spatial coordinates of the 
control points and the unknowns were simulated. Model IV 
gave slightly smaller RMS differences for the 6 unknowns and 
it was chosen for the transformation of coordinates in the 


hammer throwing analysis. 


B. Hammer Throwing Analysis 


Subjects 
The three medalists of the 1982 European Championship 


were selected as subjects for this study and their best 
throw in this competition were analysed. The mass and height 


of each subject and the distance of the analyzed throws are 


presented in Table 2. 
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TABLE 2 


MASS, HEIGHT AND THROWING DISTANCE OF THE ANALYZED ATHLETES 
Subject# Mass(kg.) Height(m.) Distance(m.) 


1 110 Ize 81.66 
2 100 94 79.44 
3 95 1330 78.66 


Data Collection 

The data were collected in Athens, Greece, in September 
1982 during the XIIIth European Athletics Championship. The 
location of the cameras can be seen in Figure 4. The cameras 
were set to an operating speed of 100 fps. (The actual frame 
rate was found to be 99 fps) An external light generator was 
used to flash a light inside each camera which was recorded 
on the side of the film. The frequency of this light was 10 
Hz. The film was Kodak Ektachrome 2537, 100 ASA. The £ stop 


was f = 5.6 and the exposure time was 1/1200 second. 
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Figure 4: Location of the Cameras 
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After the cameras were set, the calibration tree 
(Figure 3) was placed in the throwing circle and filmed. The 
set-up and the filming of the calibration tree was finished 
before the competition started. A second check was made 
after the tree had been filmed. After the filming for 
calibration, the tree was dismounted and removed from the 


field. 


Data Reduction. 

The following hardware was used in the data reduction 
process: a) A Triad Model V/R-100 mirror projector which 
magnifies the projected frame; b) A Bendix digitizing board 
with dimensions 1.215 x 0.91 m. The precision of this board 
is +0.0254 cm; c) An HP 9864A digitizer; d) An HP 9825B desk 
top micro-computer with a RAM capacity of 64 K bytes. The 
riimsprogéctorpandt &hendigit@azingoboard were*set@paratlel@to 
each other and remained in the same position until the 
digitizing of the control points and the subjects was 
completed. 

Fach film with the control points was projected and 
each point was digitized three times with respect to the 
origin of the tree. The origin from each film was marked on 
the digitizing board to be used later for the digitizing of 
the subjects. The mean values of the U, V digitized 
coordinates of each control point were stored. 

The DLT equations together with model Iv for image 


refinement were used for the calculation of the calibration 
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coefficients. The following were the two equations which 


composed model IV for image refinement: 


AUe= n@n—matee dsViex (ki mh ohrke oh eh less) 
he eP pun? ti c2octuh eP RAPE xy, C5) 
AV isu, Thee seb ioe ae fk faked ile GS ot ak aur) 
HeBaeGy? (edt +4 S2Pi xy (6) 
where: AU,AV are image refinement components of lens 
distortion and film deformation; 
a; (i=1,...,6) are coefficients of film deformation; 
U,V are the digitized coordinates of the point; 
xX = U - Us and y = V - Vs with Us & Vs the image 
coordinates of the point of symmetry; 
k, (i=1,.¢..,3) are coefficients of symmetrical lens 
distortion; 
r is the length of the vector from the point of 
symmetry to the point under consideration; 
P,, Pz are coefficients of asymmetrical lens 
distortion. 
Digitizing of the Subjects. 

Nineteen points on the body of the athlete and on the 
hammer were digitized for each subject in each frane. The 
number of frames digitized for each subject was 155 for the 
first subject, 222 frames for the second subject and 194 for 
the third subject. 

All the points were digitized with respect to the 


Origin of the calibration tree. 


Data Interpolation. 


The light which was recorded on the sides of the films 
was used to check the timing of the digitized frames from 
both films. It was found that at the beginning of each run 


the slave camera had a delay running time with respect to 
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the master camera. In order to utilize these first frames 
interpolation polynomials were used to time-match the 
frames. A first interpolation was used to calculate the 
timing of the non-matched frames. A second interpolation was 
used to calculate the time-matched digitized coordinates of 
the second film. The interpolation formula was a cubic 
Lagrangian polynomial (Gerald 1980, p.174). 
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where: P(t) is the time matched coordinate; . 
tapes, ts; ta are the times ore the currounding 


frames; yeyte 
ft1), £(2), £(3), £(4)>are the digitized coordinates 


of the same frames. 


After the interpolation, the digitized coordinates of the 
points and the calibration coefficients of each camera were 


introduced into the system of equations composed by 


equations (3), (4), (5) and (6) for each camera and the 


Spatial coordinates of the points were simulated. 
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Rotation and Translation of the Axes. 

The spatial coordinates of all points for analysis were 
measured with respect to the origin of the calibration tree. 
The center of the throwing circle at ground level was chosen 
as the origin of the reference frame for analysis of the 
throws. The coordinates of this point with respect to the 
Perdgin O of the calibration tree were X = 73.75, Y = 73.75 
and Z = -11.5 cm. The general formula for this 


transformation was the following: 
CG: = Cro ¥ Gs (8) 


where: C;,’ is the coordinate of the point of interest with 
respect to the new origin O'; 
C,o is the coordinate of the point with respect to 
the old origin O; 
Co is the coordinate of the new origin O' with 
respect to the old origin O. 


In order to present the direction of the positive Y 
axis as the direction of the throw the coordinates were 
rotated 95° about the Z-axis. The following matrix equation 


was used for the above rotation: 


x 211 4212 413 K 
Y'| = (@941 @22:423 Y (9) 
Za 231 4232 433 Z 


Boece ss X", Y', Z' are the coordinates of a point in the new 


coordinate system O'. . 
X, Y, Z are the coordinates of the point in the old 


coordinate system O. 
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OG 4 ,;SQy2;7 @ 
axis (j= 7 
the old coo 


i3 are the direction cosines of the x' 
Lede axis. (ji=-2)-and Z axisoemero porn 
rdinate system O'. 


Calculation of Error and Data Smoothing. 

The error involved in the data for the present study 
can be expressed as systematic and random error. 

The systematic error is that error, inherent in the 
coordinates of any point, which is due to the DLT method 
including the model IV for image refinement. Six points of 
the calibration tree were treated as unknowns and the RMS 
differences of the measured versus the simulated coordinates 
of these points gave an estimation of the amount of error. 
Although the error due the DLT is a random error when the 
calibration coefficients are calculated, the final value is 
a systematic error in the transformation of the coordinates 
of each point for analysis. 

The random error involved in the data was due to the 
digitizing process, that is the difficulty of detecting the 
excact position of the points of interest on the projected 
film. For the calculation of this error, a set of 7 frames 
was randomly selected and redigitized after the digitizing 
had finished. 

A second order double pass Butterworth filter was used 
to reduce the error involved in the data for analysis. The 
data were expressed in Fourier series and the power spectra 
were plotted. After investigation of the power spectra plot 


the cut-off frequency was decided to be 9 Hz. for the hammer 
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kinematics and 6 Hz. for all the other analysis. 


Analysis of Data. 


The following is a presentation of the procedures 


undertaken for the analysis of the data. 


Calculation of the Center of Mass. 


The coordinates of the center of mass of the body of 


the athletes (CMb) and the system:athlete + hammer (CMs), 


were calculated by using the following equations: 


where 


Cadaver data (Demster, 


Cb; = = [Pj ~ (Baes D; ;)qi jm; (10) 


Cs; = (M * Cb; + 7.257 * H,)/(M + 7.257) (i) 


respectively in the jth frame; 

P,;,; and D;,; are the coordinates of the proximal and 
distal end point of the ith segent; ; 

H; is the coordinate of the hammer; qi is the 
distance of the center of mass of the /th segment 
from its proximal end point in percentage of the 
botalalengthoof tthensegment; 

m, is the mass of the /th segment in percentage of 
the total body mass of the athlete; 

M is the mass of the body of the athlete; 

7.257 kg. is the mass of the hammer. The mass of the 
cable and the grip of the hammer was assumed to be © 
concentrated in the head of the hammer although it is 


in the niebourhood of 100 gr. 


: Cb; and Cs; are the coordinates of the CMb and CMs 


The quantities q; and m; were taken from Dempster's 


1955) and are presented in Table 3. 


Although these anthropometric data are widely used in 


biomechanical research, they are a source of error, since it 
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is unlikely that the body composition of Dempster's cadavers 
is the same as that of athletes and especially of throwers. 
It was decided to accept this error as it would be constant 


in all the frames and will have no serious effect on the 


results. 
TABLE 3 
ANTHROPOMETRIC DATA OF THE SEGMENTS 
From Dempster(1955) 

Segment (i) Length(q; ) Mass (m; ) 
Trunk 50.00 49.50 
Head + Neck S62 70 790 
Upper Arm 43.60 2r70 
Lower Arm + Hand 67.70 2526 
Thigh 43.30 1OV2Z0 
Lower Leg 43.30 4.65 
Foot 42.90 1.46 


(q; % of segmental length from proximal point, m; % of total 


body mass.) 


Linear Kinematics. 
The linear velocity and acceleration of the hammer, the 
CMb and the CMs were calculated numerically from the time 
displacement data. First central difference method was used 


for the above purpose, (Miller et al, 1973). 


Instantaneous Radius and Center of Curvature. 
The instantaneous radius(p) of curvature and the 
Coordinates of the center of curvature (Cx, Cy, Cz) of the 


orbit of the hammer, the CMb and the CMs were calculated by 


using the following equations: 
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v2 
SS (12) 
By — xu)? + (yy = y y)eee (2’y - zhyyeds 72 


C, = x + ———_ LENS!) 


@) 
< 
i} 
+ 


(14) 


C, = z+ —— (15) 


with : 


where : p is the radius of curvature of the orbit; 
Cx, Cy, Cz, are the coordinates of the center of 
curvature; 
xX, y, z are the coordinates of the point of interest; 
x', y', z' are the velocities of the point in the X, 
Y and Z axis; 

x", y’, z” are the accelerations of the point in the 

Same axes; 

Beis the curvature of the orbit which ts x = 1/p 


The derivation of the above equations is based on 
vector analysis and is presented in appendix C. The 


Meantitmesixhr, yigla!,ixtpex by zihwerengatculaged 


numerically from the time displacement data of the point. 


Tension of the Cable 
The radius of curvature, the linear velocity, the mass 
of the hammer and the angle of the cable with the radius of 


Curvature were utilized for the calculation of the tension 
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in the cable of the hammer. Thus, 


meu? 
Qe (16) 
p-cosé@ 


where: is the mass of the hammer; 


m 
u is the resultant velocity of the hammer; 
p is the radius of curvature; 

6 the angle of the cable with the radius p. 


Angular Kinematics 

Absolute and relative angular measurements were 
achieved through vector identities. Segments were expressed 
in terms of vectors in the three dimensional space. 

The angle formed by two segments was found by applying 
the dot product between the vectors representing the 
segments. 

Direction cosines of the vector segment were used to 
calculate the angle of the segment with the inertial axes. 

The following angles were calculated: 

Left knee angle, formed by the left thigh and the left leg; 
Hip - shoulder angle, formed by the line of the hips and the 
line of the shoulder; 

Trunk - vertical which is the angle of the trunk with the 
vertical axis; 

Radius - cable angle which is the angle between the radius 


of curvature and the cable of the hammer; 
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Trunk - arms angle which is the angle between the trunk and 
the arms (origin of the arms was taken the middle point of 


the shoulders). 


Moments and Products of Inertia. 

In each frame the moments of inertia of the system were 
eelculated about an orthogonal system of axes C(x', Y¥'., 2") 
with origin at the CMs. This system of axes was parallel to 
the space-fixed system of axes O(X, Y, Z). (See Appendix B) 

First, the moments and products of inertia of each 
segment about their center of mass was calculated for each 


frame. The following are the equations used for this 


purpose: 

e= Lee, S atel the wants Aall’ Shea ne? (17) 
Meet, wi Oi2"? tly: t@227 + 1 gi *ea27 (18) 
Mae 8 Og + 1 yp tee? + 2 aa tags (19) 
I" xy = Thaw aioel xi + @21°A22°1 yj; x M31°A32°1",;) (20) 
ie = Calta pall Mover td, pFazZ Bey h tro Take SeMesy) (241) 
eye = —(a12°a 73°]: x Pet avniaagtl’ J P9GIds2°d35-1 2 ;) C220 


mnere :risaj, I’ yi, 1°21 aré.thééprincipabsmomentsios 


inertia of the jth segment; ~ . 
Qi i, i2, @13 are the direction cosines of the 


principal axes, X (ji = 1), Y (i = 2) andOzuGis= 3) 
axis relative to the space fixed system of axes. 


The principal moments were taken from Whitsett's 


anthropometric data (Whitsett, 1963) and are presented in 
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Table 4. These data were normalized according to the mass 
and height of each athlete. The equations for the above 
normalization were taken from Dapena (1978) and are the 
following: 


For, the transverse and frontal axis, 


I = [*M*S? / MxS? (23) 


and for the longitudinal axis, 


I = [*M?%*S / M?xS (24) 


where: J, is the principal moments presented by Whitsett; 
M, S are the average mass and height of Whitsett's 
subjects M = 74.2 kg., S ='1.755 m:); 
M and S are the mass and the height of the subject at 
the present study. 


TABLE 4 
PRINCIPAL MOMENTS OF INERTIA OF THE SEGMENTS 


From Whitsett (1963) 


Segment Lice aa ies 

Trunk 1.2606 Tats oo Ossie 
Head + Neck 0.0248 0.0248 0.0168 
Upper Arm Os 02213 O02 3 0.0024 
Lower Arm + Hand 0.008:1 0.0081 0.0016 
Thigh 0. 1052 02 1052 0.0209 
Lower Leg 0.0505 0.0505 0.0050 
Foot 028036 0.0038 0.0008 


(presented measurements in kg-m?) 


(I”, principal moment about frontal axis) 
(I”, principal moment about transverse axis) 
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(I°, principal moment about longitudinal axis) 


The longitudinal axis of a segments was taken as the Z 
principal axis with origin at the center of mass of the 
segment and direction towards the proximal end point. The xX 
principal axis of the trunk was taken as principal axis 
Peraliel to the axis.of the hips and with direction towards 
the right hip. For the non-trunk segments an assumption was 
made, that is, the segment moves in a plane in two 
SuccesSive frames. The plane was defined by the proximal end 
point of the segment at the (j)th frame, the center of mass 
of the segment at the same frame and the center of mass of 
the segment at the (j+/)th frame. After the definition of 
the plane, the vector of the X principal axis was defined as 
a vector normal to the Z axis at its origin and having 
direction towards the direction of the motion. The vector of 
the Y principal axis was then defined by the cross product 
of the Z and X vectors. The method is presented in appendix 
B. 

The total moments and products of inertia of the system 


about its center of mass and for each frame, were found with 


the following equations: 


me ot Bie Bome(tyat * Cea (25) 
oH 


ly = Slity + me (Cx 3:7 + Coal (26) 
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13 

ts ‘oti. ‘BS yishees. sooutc t¢%) | (27) 
13 

ee y = ZIT xy = mE Cue (28) 
13 

Lge = ZED ei x2 m Cor Goa) (29) 
w3 

Iyz = DEI" yz: ‘hig Cy Cra l (30) 


mere: 1,, l,, 1, are the total moments and products of 
inertia of the system about the axes X, Y, 2 
m; 1S the mass of the jth segment. 
Sve sicyi) Co, are the coordinates of. the: center of 
mass of the ith segment relative to the center of 
mass of the system. 
By gg aye, CO cAtAch' seecnck's. af, the, serawetthea maments 
and products of inertia of the ith segment about the 
axes eX they ine". 
Angular momentum. 

The angular momentum of the system, was calculated with 
respect to an orthogonal system of axes passing from the 
center of mass of the system and parallel to the space-fixed 
axes. The total angular momentum of the system was found by 


the following equations: 
3 ww 
Lx = =[ Tag + L x i ] (31) 
= aM 
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where: Lx, Ly, L, are vectors representing the angular 
momentum of the system about its center of mass; 
exis lr yi 24 fare vectors representing the angular 
momentum of the jth segment about the center of mass 
of the system, considering the segment as a point 
mass; 
boent Lhyi ith eenate themvectors representing the 
angular momentum of the jth segment about its center 


of mass. 
Boeeune L',;, L',;,,~b'2; the following general equation was 
used: 

Lo =om*(r, xer2)/At (34) 


where: m is the mass of the segment; 
r,, Y2 are position vectors of the segment in two 
successive frames; 
At is the time interval between the two frames. 


me, Caw. e@and Bogasewere calculated by utilizing the 
moments and products of inertia of the segments about their 
center of mass and their angular velocity about their 
principal axes. This method is different from the one 
presented by Dapena (1978). Dapena calculated the local 
angular momentum of the non-trunk segment only about the 
transverse axis and assumed that the local angular momentum 
about the other axes was too small to be considered. The 
local angular momenta of the segments for this study were 


calculated by using the following equations: 
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where: tm iG aor en and Lie Dyer ey beat are the 
moments and products of inertia of the jth segment 
about the its center of mass; 
Wxji, Wyi, Wy; are the angular velocities of the 
segment about the space fixed axes. 


(See appendix B for details of this technique). 


After the angular momentum of the system about the X, Y and 
Z axis was found, the total angular momentum of the system 


was calculated from the equation: 


OE GO leer ae Pa a (38) 


Axis of Momentum. 
Axis of momentum of the system for each frame was 
defined as the axis which passes through the center of mass 
of the system and is coincident with the angular momentum 


vector. The direction cosines of the angular momentum are: 


Gees) Gimeno ly? ue Ee (39) 
Pee (eet Ly? + Le ee (40) 
6) ct ON (a a se ann ice eb (41) 


The coordinates of the point of intersection of this axis 


with the ground are: 
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A= (65%, Z1°B1)/B3 (42) 
Yoo OBig = Vi a3 Z1°B2)/B3 (43) 
ZO (44) 


where: X,, Y,, Z, are the coordinates of the center of mass. 


Instantaneous angular velocity. 

The instantaneous angular velocity of the system was 
defined as the angular velocity of the system about the 
System of axes passing through center of mass and parallel 
to the space-fixed system of axes. The following matrix 


equation was used for this purpose: 


ie Py Ixy I x2 2, 
Ly ea ig a cael Qy (45) 
L, -Ix. = viz I, Q, 


Mnere: 9,, Q,, 2, are the angular velocities about the x, ¥, 


Z axes; 
Pei) vl stick yyy. x cae vines ee OMelC’ and products of 


inertia calculated with equations (17) to (22); 
L,, Ly, L, are the angular momenta of the system 
calculated with equations (31) to (33). 


The column matrix of the angular velocity in the equation 
(45) was found by multiplying the inverse of the moments of 
inertia matrix by the column matrix of the angular momenta. 
The magnitude of the resultant instantaneous angular 


velocity of the system was found by the equation: 
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v= -( Q,2 + Q, ? = Q,72)172 (46) 


External Torques of rotation. 
The external torque of the rotation was defined as the 


rate of change of the angular momentum. Thus, 


 =SAL0Y Wt (47) 


Forces of Translation. 


The external forces of translation applied to the 


System were found by using the general formula: 


(48) 
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ll 
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where: m is the mass of the system; 


a is the acceleration of the’CMs. 


Kinetic Energy. 


The kinetic energy of the system was found from the sum 


of the kinetic energy of rotation plus the kinetic energy of 


Mnanstation. Thus, 


Tis Te -+ Tr (49) 


where: Tt is the kinetic energy of translation and Tr is the 


kinetic energy of rotation given as follows: 
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Tt = ——M-v? (50) 
2 
1 
Be LO. + Ieee + ie, 7 
2 
TARPS LOA BNA KiOFos Os DE, shyo gy) CS) 


where: M is the mass of the system; 
u is the velocity of the center of mass; 
Pal. he and thn, Teoy, iyo are: MOMeENtS wancuproaucrs 
of inertia of the system; 
Wy, Wy, W, are the angular velocities. 


All the programs for the data reduction and data 
analysis were written in HPL computer language by the author 
and the HP 9825B desk-top computer was used. The motion was 
Givided into turns and each turn into single and double 
Support phase. The results of each phase, each turn, and, 


each subject were compared to one another. 
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IV. RESULTS AND DISCUSSION 
The sequence of presentation of the results found in this 
study 1s as follows: Results of the DLT tests; Error 
analysis; Kinematics and kinetics of the hammer; Kinematics 
and kinetics of the body Center of Mass; Segmental 


kinematics; Kinematics and kinetics of the syste m. 


A. Results of the DLT Tests. 

The DLT method together with six mathematical models 
for image refinement (Karara and Abdel Aziz, 1974) and three 
different calibration trees were tested. Criterion for these 
tests were the Root Mean Square (RMS) differences of 
measured versus simulated coordinates of the control points, 
or of control points which were treated as unknowns. 

The results of tests on the first tree which was a 2-D 
frame filmed in different angular positions were 
unsatisfactory in that there was considerable deformation of 
its shape in rotation caused by its weight. 

The results of tests on the second tree which was a 3-D 
frame are presented in Table 5. 

TABLE 5 


RMS DIFFERENCES OF THE CP FOR THE SECOND TREE 


Model X RMS Y RMS Z RMS MEAN 
No Ref. iieAS 1.0214 0.8458 PSO 162 
Model II 1% 1231 1, 0998 OO he 1.0147 
Model III ime 109 1.0944 0,8929 t, 009 T 
Model IV 141329 SPM toe ges O. 7614 1.0245 
Model V et ti5 1.4940 0.7676 1.1244 
Model VI 1.3235 129989 ORS BUS 123030 


(presented measurements in centimeters) 
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The results of the third tree which was geometric in 
Shape and was used for the data collection of hammer 
throwing, are presented in Table 6. In this Table the first 
column contains the RMS of all the control points, the 
second column contains the RMS of 27 control points from the 
test when 6 points were treated as unknowns and the third 
column contains the results of the 6 points treated as 


unknowns. 


TABLE 6 


RMS DIFFERENCES OF THE CP FOR THE DATA COLLECTION TREE 


Model 33.eGP 27 CP 6 Unknowns 
No Ref. X 0022 Oe29 0.84 
Y O20 0.14 C239 
Z 0.24 0.19 O85 
Model II xX 0. 36 0.42 0.42 
Ne 0.33 C236 O55 
Z O25 Ons te OF61 
Model III X 0.24 O230 0.46 
Y 0.19 ON6 O59 
Z O8Z25 0.19 0.65 
Model IV X 0.24 ORzsS O27 
Y ORG Re 0.14 0.60 
Z O225 Oise O72 
Model V yA 0.29 0.29 Os30 
Y Oe G O218 0.74 
Z 0.25 0.20 0.86 
Model VI X 0.48 O25 1.46 
Y @. 36 0.14 0.76 
Z O. 36 0. #S oa 


(presented measurements in centimeters) 
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The digitized, measured and simulated coordinates of the CP 
of the last tree, as well as the calibration coefficients of 
the DLT with model IV, are presented in Appendix D. The RMS 
error presented in Table 6, was the error of the DLT method 
in the present study. In a similar study of the DLT method, 
Shapiro (1978) found an error of +0.40 cm. in the X-axis, 
20°40 7cm. in the Yramishand»1+0.50;°cm. win the Z-axis.» Dapena 
et al (1981), found a RMS(X) = 1.50 cm., RMS(Y) = 1.32 cm. 
and RMS(Z) = 0.60 cm.. With the second tree, Model III gave 
Slightly better results than the other models, while with 
the data collection tree, Model IV gave the best results. 
Another test on the data collection tree was the 

Simulation of the coordinates of the location of the cameras 
with respect to the origin of the tree. These coordinates 
and the simulated distance of the cameras (£4, and £2) from 


the origin are presented in Table 7. 


TABLE 7 
SIMULATED COORDINATES AND DISTANCE OF CAMERAS 


FROM ORIGIN WITH DIFFERENT MODELS 


Model Camera #1 Camera #2 
X Y Z £3 X Y Z L> 
No Ref. 358758542100 n70c386x13 37 «64424598 0292 45.16 


Model II 35.86 -5.22 0.70 36.24 37.12 25,0400. 92 425,26 
Model III 35.95 -5.23 0.70 36.33 37268425. 01m Ors Zee 5ues 
Model IV B59 -5e 22. 0.10 36.22 37/60 220002 Un oe, 20.20 
Model V 35.87 -5.18 0.70 36.25 37 6582550000592 45.20 
Model VI 86040) —5.26.0,7 1) 36.70 30022) 25.4000. 72eeo,06 


(presented measurements in meters) 
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The X, Y, Z coordinates of the location of the two 
cameras with respect to the origin were not measured 
directly in the filming process. Only the distances of the 
cameras to the center of the throwing circle were measured, 
and they were found to be £, ~ 36.20 m. for camera #1 and 
ke, =~ 45.10 m. for camera #2. 

The results from the third tree were superior to the 
ones of the preliminary tests because in the preliminary 
tests the CP were randomly distributed in space and it is 
possible that they included higher measuring error. 

In order to test the stability of the DLT for outside 
areas of the calibration tree, 6 CP located in this area 
were treated as unknowns (points 28 to 33 in Figure 3). The 
RMS error of the measured versus the simulated coordinates 
of these points are presented in Table 6, 3rd column. These 
results indicated that, use of a well constructed tree and 
Karara's models for image refinement enabled reliable data 
woebe Collected. usingsa’ calibration’ tree: smalilerothan the 
area of action. These results were expected, since the 
control points were used for the estimation of coefficients 
that indicate the orientation and location of the cameras in 
the filming process and the orientation of the secondary 
image in the digitizing process, rather than for the 
estimation of ratio parameters of control and uknown points. 
After these coefficients were computed, the spatial 
coordinates of a point under consideration were calculated 


based on the digitized coordinates of the point and the 
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calibration coefficients. Utilization of the mathematical 
models for image refinement corrected the transformed 
coordinates of a point if they were affected by film 
deformation and/or lens distortion error. These models were 
of significant value when a point under observation was 
located in the outer area of the image in film from at least 
one camera. Points near the outer area of the images in 
films from both cameras were expected to include larger 
error than other points. These areas however, were not used 
for analysis because a larger area than that taken up by the 


subject was filmed. 


B. Error Analysis 

The error involved in the analysed data of hammer 
throwing was classified as systematic error or random error. 

Systematic error was the error inherent in the use of 
the DLT method. The main sources of this error were: a) the 
least square method used to derive the calibration 
coefficients: b) the limited precision of measurement of the 
CP's spatial coordinates and c) the digitizing error of 
these points. Another source of systematic error was the 
anthropometric data used for the length, mass and moments of 
inertia of the segments. The criterion for the estimation of 
the systematic error was the RMS difference of the measured 
versus simulated coordinates of the CP. This error was found 
to be RMS(X) = 0.24 cm., RMS(Y) = 0.19 cm. and 


BMoeZ) a= 0.25 cm... 
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Random error was defined as the error mainly created by 
the researcher's inability to detect the exact location of 
points in the digitizing process. Other sources of this 
error were vibration of the cameras in the filming process 
and vibration of the projector and digitizer in the 
digitizing process. This error was estimated by redigitizing 
seven randomly selected frames. The RMS differences of all 
the points from the two digitizings were found to be 
RMS(X) = 1.73 cm., RMS(Y) = 0.64 cm.,; RMS(Z).= 0.62 cm..:The 
differences were higher for poorly detected points and 
smaller for well detected points. The absolute differences 
of the two digitizings for the right hip (poorly detected 
point) and the hammer (well detected point) for the seven 


frames are presented in Table 8. 


TABLE 8 
ABSOLUTE DIFFERENCES OF TWO DIGITIZINGS 


FOR THE RIGHT HIP AND THE HAMMER 


RIGHT SHIP HAMMER 
Xx vy Z x y Z 
FRAME #1 3.41 PRS aS ese ew: OAC Gl aUrsc 


FRAME #2 3.34 Cake 10466 Oaze Oie7 HKhO228 
FRAME #3 0.66 Gua]. 18 Uo O77 08-8 05 40 
FRAME #4 1.04 OF46 pia86 0501 0 200530 £0 
FRAME #5 Paes 74 2305) 0.26 0.20 Qe OameU aon 
FRAME #6 2.56 beO 5.9 [0:531.2 0.76 QO. Doe. 22 
FRAME #7 0.48 1.69 0.48 1.84 OnO7e > OnZs 


(presented measurements in centimeters) 
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It should be noted here that the X, Y, Z-axes presented in 
the above tests were the axes before the rotation, therefore 
the X-axis of the tests is the axis of the direction of the 
throw in the analyzed data. 

Digital filters were used for the reduction of error in 
the data. The cut off frequency used was Fc = 9 Hz. for the 
hammer kinematics and Fc = 6 Hz. for all other parameters. 
This decision was based on study of the power spectra plots 
of the signal. Figures 5 and 6, are two examples of these 
plots for the X-axis of the hammer and the right hip points. 
An example of raw versus filtered data can be seen in Figure 
7 which shows the acceleration of the hammer. This 
acceleration was calculated with first central differences 
from raw data and data filtered with cut off frequency Fc = 


9 Hz. 


C. Kinematics and Kinetics of the Hammer. 

The analyzed throws were (a) a three turn throw of 
mooem. (Subject 1); (b): a four turn throw of we.44 Mm. 
(Subject 2); and (c) a four turn throw of 78.66 m. (Subject 
3). Each turn was divided into single support phase and 
double support phase. The beginning of the single support 
phase (BSS) and the beginning of the double support phase 
(BDS) for each turn were used as reference instances for the 
analysis and presentation of data. Other reference instances 


for each turn were the times when the hammer reached the 


highest (MAX-Z) and the lowest (MIN-Z) points of its orbit. 
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Acceleration of the Hammer, Raw and Filtered Data, Subject 1 


Figure 7 
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The analysis of each throw started at the beginning of the 
first turn. No analysis was made of the preliminary Swings 
of the hammer. 

The height and angle of release were found to be 1.27 
Bmmences0"/? for Subject-1, 1.26 m?,and 37.3° for PSubject 2 
Spemiecs M. and 37.6° for Subject 3. The angles of release 
for all the subjects were less than the calculated optimal 
angle of release which for an 80 m. throw was 6 = 44.3°., 

The maximum velocity of the hammer was found to be 
Bees measmfor Subject’ 1,930.72 m/sifor Subject i24and 31.02 
m/s for Subject 3. The maximum velocity of the hammer 
occurred 10 to 20 milliseconds before release for all 
Subjects. 

The timing of each throw and reference instance are 
presented in Table 9. In each turn the hammer reached the 
MAX-Z point in the single support phase and the MIN-Z in the 
double support phase. The time of the single support phase 
was shorter than the time of the double support phase for 
all the subjects and all the phases. Subject 3 had the 
fastest turn time for each turn. The same subject spent less 
time than the others in driving the hammer upward during the 
double support phase 22%, 27%, 24%, 38% of the double 
Support time for each turn. The eqivalent percentages for 


Subject 1 were 29%, 30%, 46% and for Subject 2, 31%, 32%, 


22k, 39%. 
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TABLE 9 


TIMING OF EACH REFERENCE INSTANCE 


Point subject “Subject 2° Subjects 
MAX-Z Sars 0.29 OR aks! 
BDS oe 0.41 0.24 
MIN-Z sos 0.68 0.49 
BSS i 0.80 OFo6 
MAX-Z 0.23 98 O75 
BDS UR oe 107 Orr 
MIN-Z Oro] Wee 2 O96 
BSS 0:59 e229 ibe is) 
MAX-Z Cero 1.45 1320 
BDS Do) Oe Paz 
MIN-Z 04 9:7 ital sy 1.41 
Doo 1.04 ieer3 1.46 
MAX-Z bean e386 ei6.3 
BDS hee RS Te) 1x67 
MIN-Z 1,39 2.09 Peeks: 
RELEASE 1.51 218 pe 


(presented measurements in seconds) 


The spatial coordinates of the hammer for the different 
Subjects showed that the hammer traveled in a spiral mode 
starting at an acute angle to the horizontal and finishing 
mith thearnelease angles, inthe’ first turn, Subjecti2 rotated 
the hammer in a plane closer to the horizontal i.e. at a 
more acute angle than the other subjects and gradually 
changed this angle to the release angle. Figures 8, 9 and 10 
present the trajectory of the hammer in the different planes 
FOL subject 3. 

The radius of curvature of the hammer for each of the 
three subjects can be seen in Figures 11, 12 and 13. The 


length of the radii oscillated from phase to phase and 
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within the same phase. There was a Gradual decrease of the 
amplitude of oscillation from turn to turn. For each turn 
the length of the radii became maximum at the beginning of 
the single support phase, gradually decreased during the 
Single support phase and began increasing again in the 
double support phase. In the last turn, and after the radii 
had fallen to their minimum length, the amplitude of 
oscillation became very small for Subject 1 and 2. Subject 3 
increased the radius after the BDS point and decreased it to 
a minimum length when the hammer was at the MIN-Z point. In 
comparison with other subjects, Subject 1 presented the 
smoothest oscillations with the highest amplitude, while 
Subject 2 presented the noisiest oscillations and Subject 3 
the smallest amplitude. In the last double support phase, 
after the hammer passed from the MIN-Z, the radii's length 
increased and approached infinity at the release point for 
Subjects 1 and 3, while Subject 2-released the hammer before 
Meestarted traveling in a straight ‘bine (1.¢e. with low 
length of radius). 

The angle between the radius of curvature of the hammer 
and the cable for all the subjects ranged from 0° to 13° 
Figure 14 represents the change of this angle forrSubject 7, 
while Subjects 2 and 3 showed similar curves. This angle was 
greater in the first turn and decreased gradually until the 
Second from the last turn, then increased again. During the 
single support phases the angle started with a small value 


and increased until the middle point of the phase. Before 


. uex 
it : : i ree 
3 | Be La aN 
os a 
eit Io sas6r seb Ltueete 6 new e4 


regs (sas 10% eeu as ara7 air 
5< online gad aad JR Met ee snieaed vas oi 
od’: pring Beeee teed yf ap BRIS Sear 
Sa2 Wh er Sica est os iii fine bully 
‘Tas? adv 29726 "DG6 base ted reel ‘oat ol: .ssad@ 2 
su St tometer sete) an sigdg o4ad 


f suatdor 8 bae, Taeetee 102 (5 begs ‘yiay, page: 
dakdanee pie) 2atcd AoE $52 yeitte austew. . 
gare, Sa ge3 $h bew “ister ent nadie ot ; 
ech b4yngesag)! Aestaur . 3328 RRe  Tsar0 bgt 
Ssictiw pry beflomte tc etn 4 ais iv maakanttiaee # 


i ssetdvd Sadcensiseiiitac sesneggn eR tty 


a TIAN ‘ 


ame ticeqye, eldued s20l sae tus gn 
daptie Fe’ EF ppt air, co- M) otd mips 


- ~~ = A 


sr) 90. 38 Sduoither sae 7s4 sist ha ope yi 
610is¢  :9mter oan? Ostaeiaa Sits ate hid t fos a 


vated .2f> Bo ‘pyganaa>: ec as-cnauee 


iris: 


——— ¢ 7 | 7.4 
“t Seetiive sot Sia wins 20° aes. ; 


anata «9ouigQ sii 


72 


| 30efqns ‘uewweH 40 91qeQ ULM euNn}eAUND 4O SnNLpey jo oe Buy 1 eunbBig 


(99s) GWIL ; 
[3°0 6S°0 ZE*0 


tT 


oT 


6T 


be 


ba 
asva1a¥ Si , (Dep) 


c2za0ow 


Ba 


i 


THe aBMIc 


Toet@ue , TS98hh 


i 


a 
Swe 


~ 
4 
A « 
i 
i s 
a BN me er me 
, 


i rs 7 
: 7 | 
g SEA 


$ 
i 


= 


eS ee ee eee Pe 


a = 


2 o 
“23 | : S 


i362! ~y 


30 ES (Hai osaecss 


i 


2 


7 
‘ 


q: 
. 


“vsuanele ajavod.20 aviewigalesas 


wo to eyiben Jo sige i0f emagra 


73 


the hammer arrived at its MAX-Z point the angle started 
decreasing again arriving at a minimum at the end of the 
Single support phase. In the double support phases the angle 
decreased from the beginning to the middle point of the 
phase and then increased again. In the last two turns and 
during the transition from double to single support and from 
Single to double support the angle decreased. In the last 
turn the angle started increasing after the MIN-Z point 
arrived to a maximum angle then decreased until = 40 
milliseconds before the release and started increasing again 
until the release point. 

Figures 15, 16 and 17 represent the velocity curves of 
the hammer. The velocity increased in the double support 
phases and decreased in the single support phases. There was 
@ gradual increase of velocity from turn to turn. The 
contribution of each turn to the maximum velocity can be 


seen in Table 10. 


TABLE 10 
CONTRIBUTION OF EACH TURN TO THE MAXIMUM 


VELOCITY AND ACCELERATION OF THE HAMMER 


TURN VELOCITY ACCELERATION 
Sat SoZ Sree: Soi Sa2 S.2 
First es 69 76 ee 56 61 
Second te 80 83 fie yey 80 
Third 86 86 89 86 80 89 
Fourth 100 100 100 100 100 100 


(presented measurements in percentage of maximum) 
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There was a local minimum velocity at the BDS POINne an 
each turn, except in the second turn of Subjects 1 and 2 who 
achieved the local minimum velocity before that point 
(Figures 15 and 16). In each turn the hammer reached its 
local maximum velocity at the end of the double support 
phase and beginning of the single support phase. In the 
first single support phase the velocity curves showed great 
oscillations, probably caused by the athletes' effort to 
achieve a high velocity immediately. It was the velocity in 
the Y and Z axes which primarily contributed to the local 
maximum velocity of the hammer in each turn. 

In the single support phases for all the subjects and 
for all but the first turn, the decrease in velocity was 
calculated in percentage of the maximum velocity. This 
decrease was found to be: 6%, 10% in each turn for Subject 
eC ope Gm in—each=turn fLor-Subject-2-and-8%, 77%, 12% in 
Bacuecurn for Subject 3. 

The accelerations of the hammer for all subjects are 
presented in Figures 18, 19 and.20. The hammer accelerated 
in the double support phase and there was a gradual increase 
from turn to turn for each throw. The highest acceleration 
for each turn occurred at the time when the hammer reached 
the MIN-Z point of its orbit. The maximum acceleration was 
found to be 478 m/s? for Subject 1, 505 m/s’? for Subject 2 
and 492 m/s? for Subject 3. This global maximum occurred 60 
milliseconds after the local MIN-2 point of the hammer and 


60 milliseconds before the release point for Subject 1. The 
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eqivalent times were 50 and 40 milliseconds for Subject 2 
and 0 and 80 milliseconds for Subject 3. Subject 1 kept a 
high acceleration at the global maximum for a longer time 
than the other subjects, while Subject 3, after reaching the 
global maximum presented a second maximum of lower value 
approximately 50 milliseconds after the first peak. The 
contribution of each turn to the maximum acceleration for 
all subjects can be seen in Table 10. It was mainly the 
acceleration in the Y and Z-axes that contributed to the 
peaks of the resultant acceleration (Figure 21). In the 
Single support phases the acceleration increased slightly, 
in the second turn for Subject 1, in the second and third 
murnetor Subject™2.andg, in the third turn for Subject 3. ein 
all the other single support phases the acceleration either 
remained the same or decreased. The accelerations presented 
different oscillations for the different subjects. Subject 
1, increased the acceleration in the second single support 
Dnase, while, in the last single support the acceleration 
decreased until approximately 30 milliseconds before the 
hammer's MAX-Z. Then, it increased until the MAX-Z and when 
the hammer started traveling downwards to the BDS point the 
acceleration decreased again (Figure 18). Subject 2, in all 
the single support phases increased the acceleration from 
the BSS point until approximately 20 milliseconds before the 
MAX-Z. Then, he decreased the acceleration untid the MAX-Z 
and increased it from there to the BDS point (Figure 19). 


Subject 3, in the second single support phase showed 
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opposite oscillations to the ones of Subject 2. In the last 
two turns he decreased the acceleration until approximately 
10 milliseconds after the MAX-Z and increased from there to 
BRE"BDS point (Figure \20). 

Figures 22, 23 and 24 represent the tension in the 
cable of the hammer calculated using equation (16). The 
tension was maximum when the hammer reached the lowest point 
of its orbit. This maximum tension was found to be 3730 N 
BoresuUbject 1, 3602 N for Subject 2 and 3837 N for Subject 
3. The contribution of each turn to the maximum tension was 
momar tOube 62%,.79% and ‘100% for Subject 1) 50%) 72%, Giz 
meri, for Subject,2 and, 55%, 72%, 82% and, 100% for 
Subject 3. For each turn the local maximum tension occurred 
in the double support phase. In the single support phase the 
tension decreased until the point where the hammer reached 
its MAX-Z and from that point the tension started increasing 
again. Only Subject 1 achieved an absolute increase of 
tension in the single support phase and only for the second 
e Og oe 

The impulses applied to the hammer in different phases 
and for the different subjects are presented in Table 11. In 
the last double support phase Subject 1, applied a greater 
impulse than did the other subjects. This reflected the 


application of force over a longer period of time. 
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TABLE 11 


IMPULSE APPLIED UPON THE HAMMER FOR EACH PHASE 


Phase Subjected Subject 2 Subject 3 
SS et eter SOC mL. | 243.47 
DS | Be ae a eM We Ty 524.42 
SRS) 419.01 439.807 454.73 
DS 532.96 £92533 O5'2..3/14 
So 484,44 SLecos 497.06 
DS Cae Siew! tl Ose 627.58 
Ss 556'.7S Soto pis) 433.49 
DS 744.32 TSOeS2 L205 Ge 


(presented measurements in Nxs) 


Figures 25, 26 and 27 represent the moments of inertia 
of the hammer about the orthogonal system of axes parallel 
to the inertial frame and passing through the center of mass 
of the system (CMs). In the first turn Subjects 2 and 3 
achieved larger moment of inertia about the Z-axis than did 
Subject 1. This was because they rotated the hammer in a 
plane which was closer to the horizontal than to the 
vertical plane, with Subject 2 presenting an angle with the 


char 


Im 


horizontal 6 
Pagures 28yeco™ and -30. represent the distance) of the 
hammer to the CMs point. This distance increased in the 
Single support phase and decreased in the double support 
phase. In the double support phase athletes pulled the 
hammer toward their CMs and in the single support phase the 
hammer traveled away from the CMs. In the first double 
support phase, the distance increased to the point where the 


hammer reached its MIN-Z, then it decreased again until the 
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Figure 25: Moments of Inertia of the Hammer, Subject 1 
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f the Hammer, Subject 2 
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Moments of Inert 


Figure 26 
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Figure 27 
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Hammer - CMs Distance, Subject 2 


Figure 29 
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Hammer - CMs Distance, Subject 3 


Figure 30 
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BSS point. This happened with all the subjects, but was most 
evident in Subject 1. 

Figures 31, 32 and 33.nepresent the angular momentum of 
the hammer with respect to the CMs. The angular momentum 
increased in the single support phase and decreased in the 
double support phase. There was a gradual increase of the 
angular momentum from turn to turn. The maximum angular 
momentum was found to be: 327.95 Kg*m?/s for Subject 1, 
21.44 Kg*m?/s for Subject 2 and, 319.81 Kg*m?/s ‘for Subject 
3. In the first single support phase, Subject 1 reached and 
mantained until the beginning of the second single support 
phase, a local maximum angular momentum which was 83% of the 
global maximu m. Subjects 2 and 3, reached 65% and 76% of 
the global maximum respectively and continued to increase 
the angular momentum until the end of the turn. The 
contribution of the last three turns to the maximum angular 
momentum was respectively: 83%, 94% and 100% for Subject 1, 
Bes, 93% and 100% for Subject 2 and, 88%, 96% and\ 100% for 
Subject 3. Subject 1, achieved in the end of the first turn 
what the others achieved in the end of the second turn in 
terms of angular momentum of the hammer. These results, 
compared with the radii of curvature of the hammer (Figures 
11, 12, 13) indicate that the angular velocity of the hammer 
increased in the single support phase when its radius 
decreased in length. 

Figures 34, 35 and 36 represent the angles between the 


angular momentum vector of the hammer and the X, Y and 


we 


quan 25H 


x 
: : 
fd re 


Lage ciao aad en? Lis ithe bensasqart eidt .20: . 


{ c 7 : onal 
} insmom, Isleping) ate Aaene2Qet fs bis we 12 eek 
r , ~~ “ r. ie , rs as” I od er 
anmenreanot 2S Lu0pne hae Se pee oY S28GRS2 iia 


a 
apsstpeb! bah sesdg 77 ON sigote af ant 


Janata? Lewoseiag & any vied? conask +1 Oa a 
Yelvera’ mb ai idett Sey smut 2 al aot “au JOR 3 
nepabe Ata AS ES cay Se ‘aid 62 Sned? ow aE 
; a\ free Nase e” bags jsetdue +91 ‘\" wept | 
aid eat! SasQaee vt hcg sage etente seme8 a 
. ne —) hitoaee “er ae ont ankged 4.3 tisne * 
sw: Core ML rice ‘sti Ruins & Sagan tasol a 
we 425 Bad sks 3 *s -Siee S -aeoer“ue —_ 
a pak Sapeidtib.s | wore wha ae fe 
. it. deur shi io tee any hein mu dramor ratua 
‘upie mim icon add en? ogee init afd Yo nobae 
aL i Joh | OG .OFe, ae aie piofawirsaqe>s ea | 
| 16 hae cea dene ak jo} XOGr bra 
$2969 sitio inerens | » baveides .' sandud, f ; 
i { 
y4i Goose att to ots ota ab bovetdos’ sr8dio-# 
siivass eget aeons asta ke Wy toomom satugte lo 8 
eStupi a) omih\ aaa ba 234 tavaae 20 Lites old ast, B 


1 canner! 


> 


aay « 


mine Git 2a Lons. > ‘ane ae ban ct at tA 


re 


Str 6 V tiaolev ielueae its 1603 . etentbas ery . 


sutbey att trate Saag sanaque aan. ods ~~ wig 
f . ispaod, at 


Sy 


bas YR saz be : al ae sOt90e av snsivom 


i aN 7 : : 
‘ L ao ’ ) . a 9 


2D) 


| yOefqns ‘ueuweH Sy} JO WN}UEWOW ueLNBuy :1¢ eunBLy 


LuYOddNS JISNIS 3O ONINNID3E=SSa ‘LYOddNS 37gnod 4O ONINNID3ag=Sag 


(99s) AWId 
BEE GZ" T pO°T Teo 6S°O 


0S 


OOT 


OST 


002 


a 


OS¢é 


Ke 


00t 


(S/WyWyb5y) 
aS WaTda sda Sia sdd SS. Sich Sdd 


aw 


te. 


| 
. 2 
i * 
| 


4 ams so arn Ad 
i 

' 
2 


a hoomes 


[Oe ie ei 


Say 


a , 7 & Seg 
: oS a BS Ch eS 


vine - a err _ 


é 
] = t 


» . : 
al is 7 
&., a ae SS zy a os aa Dh bit 
ba = ‘= 
ee ee 
1. . a ae 


96 


@ yOefqns ‘yuouwey ey} 4O WN}USWOW ueLNBuy :zE eunbi4 


LuYOddNS JTIONIS 4O ONINNIDAG=SS@ ‘LYOddNS 31gnod 4O ONINNID3g=sag 


(99S) AWIL 
SLeZ G6°T CL St CGE 6ccm EOE 


OOT 


O2T 


OVT 


09T 


O8T 


002 


072 


ZoztuwazkbH oe 


0v2 


O 


092 


a 


082 


4SWa Td S dd SSq | sdd SSd sdd S32 


sdd 


oO y 

eX. : er oan Ps Sake 
7 neat 
an Ps - » PCat ‘ oF 
ny , i ae 

{ Bet ’ . . 

: ) 1 ‘ j F i ° 

Ve 
it 2 ve x) j j 
7 
7 ~— 
us > 
‘ 
el ee TI eee Necro eee ih a arm ang tip—caneaaata: 
A 
ie eee ye et dn eae 
; ] PP PEIN ids eae, tp ay x 
, ; i 
4A 
rf . - 
q ‘ Pig ~ 
» ae j he, es 6 
7 — = - ~~ nt oe (ee Pao er 


OF, 


€ }oefqns ‘ueuwey ey} Jo uN} USWOW ve_nBuy :€€ eunbiy 


1YOddNS JFISNIS 4O ONINNISSG=SSE8 ‘LYOddNS 31gNoG 4O ONINNIDAg=Sag 


(99s) GWI OL 
T6°T LO°T 9b°T Gort Gost 6L°0 


Ov 


08 


OT 


09T 


SOE HW ZHD EZ 


002 


a oO 


OVC 


cc 


082 


(Ss /W,yW, by) 
aS Va Tau 


te 


eeOy i 


‘ce 


] in _ 


a) _ 
Bet 
a 


> 


ae. 


‘| 


e 


x oa aT 


Sn ennai. 


TIS, See at 
-_ 
‘a _ 


98 


L 30efqns ‘yo}SeA wWN}UEWOW ue_NBuy s , yowwey 04} 40 Sal Buy UuolL}OeuLQq ‘ve eunbiy 


(Z)SIXV-Z “(O)SIXV-A ‘(X)SIXV-X 
LYOddNS JTIONIS dO ONINNIDSE=SSE ‘LYOddNS 37Esnod 3O ONINNIDg’g=Sag 


(99S) GWIL 
yates yO: T8°O 6S°0 ZERO 


Oc 
Pas PR. Gy 


L Poeecces eatstes 7 
Z Zz “q 


¥ co cedathel tl PP Z 


, ey A en Ae tee 


ee 
Pecee © 000 0008 000 05 5590090? 


OV 
Zag 


09 


1 O8 
yk 


022800, |X x xX x 


Dera or 


x Xo." RICO OC 


O 


4S Va Tdu 


OOM © 
O 


reper een O 
e) pono n cnet 


Gee 6 
XK aererertessetses, 


: “Took ee Ss GE 3 
°° Sato cobs ck anttossaces OOT 


OcT 
O 


(5eap) 


aq_2z2orn 


=? ¢- 


ny 
5 


ii je 


Hmeci mcr F 


av 


of 

ine 

Wy 

sf | 

uh? 

* 

‘ Ae 

As 

~ a 
im) 
err 


bo a 


ith 22 oe 


: 
~ 
a 
© 


A 
=" keh 
i. o 
p< 
Me 


=F 


SP SeRSST HAIG “OF DON 


“ner 


og 


f 30efqns ‘yO}SeA WNyUEWOW UeLNBuy S ,dowweH 84} $0 se,Buy uoljoeutgq :GE aeunBL4 


(Z)SIXV-Z ‘(O)SIXV-A *(X)SIXV-X 
LYOddNS JISNIS 4O ONINNIDSE=SSg “LYOddNS 37Egnod JO ONINNID3S8=Sag 
(99S) AWI dL 


SB IEC SC. | GL I C.aeh Cceet BOAT 


08°0 


OT 


O€ 


OS 


OL 


06 


OTT 


O€T 


(fap) 
aS VaT aa 


sdqd SSa sdqd ssa sdqqd Ssqa sdqd 


ae 2own 


100 


€ 30efqns ‘4O}SeA WN}UEWOW ueLNBUuy s ,ueuWeH ay} Jo se[ Buy UOL}O8ULq :9€ eunBLy 


(Z)SIXV-Z “(O)SIXV-A ‘(X)SIXV-X 
1Y¥OddNS JTIONIS 40 ONINNISSE=SS8 ‘LYOddNS 37gsnod JO ONINNID3a=sag 
(99s) AWIL 


G Opel 6L°0 


OT 


Of 


0S 


OL 


66 


OTT 


CET 


(bap) 
aS Wa TSH 


Sdd 


“e{20H 8 


I 
tina 


L 
oo | 
oe! a om ’ 
1 { | 
‘ _ 4 4 
‘ | a a f 1 Aw aght ¢ 
ea i Set 20: 
ee | | | 
/ : } 
| ; | - | 
| | 0 aE 
b 4 
: a Seta areata ieee WA — a as 
a. ee ee A a a ee oe ree ——i nee meena ee es 
: J 
i ae pre c » BYLO 7a. 4 Rc 
m “ ' es 5 we = z ; 
TFGe ae 3 sate 3c) meTieinseeecs. .TeOwve s2eu00- 30 PA etigeees = 
- '.¥ i i i Aft= iy JS TX A > - 
— : 
S° E teetde .Apicev mutasat ~sT up-c ‘sembi ott Jo celgnk qernaaan emeaianpem 


101 


Z-axis. These figures show the angular position of the plane 
of the hammer in space, which will be perpendicular to the 
vector. The angle of this plane with the vertical at the 
beginning of each throw was 6 = 35° for Subject 1, 6 = 9° 
Boreoubject 2 and @ = 25° for Subject 3. This angle 
increased during the single support phases and decreased 


during the double support phases. 


D. Kinematics and Kinetics of the Body. 

Selected trajectories of the center of mass of the body 
(CMb) in different planes are presented in Figures 37, 38, 
eSJeand 40. In the Y-Z plane, Subjects'2: and 3 had similag 
trajectories of CMb. In the X-Y and X-Z planes all three 
Subjects presented similar trajectories. 

The maximum vertical displacement of the CMb in each 
@ucn was less than 0.12 m. for all subjects. In the Z-axis 
the CMb moved in the opposite direction to the hammer except 
during the last turn when the maximum height of the CMb 
occurred at the release point same as the maximum height of 
hammer. 

Figures 41, 42 and 43 represent the CMb radii of 
curvature. These figures indicate that the highest rotation 
of the CMb occurred at the end of the double support and 
beginning of the single support phase. The radii of the CMb 
became smallest at the same time as the radii of the hammer 
became greatest. Subjects 2 and 3 presented noisier 


oscillations of radius than Subject 1. 
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Figures 44 and 45 illustrate the resultant velocities 
of the CMb for Subjects 1 and 2. The velocity of CMb for 
Subject 3 was similar to the one of Subject 1. Subjects 1 
and 3 presented a gradual increase of CMb velocity from turn 
to turn with the global maximum in the last turn. Subject 2 
achieved the global maximum velocity in the second from the 
Hast turit=The maximum. velocity for each turn. occurred at 
the beginning of the double support phase and remained 
relatively constant until the middle of this phase. After 
this point there was a rapid decrease in the velocity and it 
was least during the transition from double support to 
Single support phase. At the same time the velocity of the 
hammer reached its local maximum for each turn. It was the 
velocity component of CMb along the Y-axis which mainly 
contributed to the maximu m. This, together with the Z 
component (Figure 46) indicated that in the first half of 
the double support phase the direction of the CMb's velocity 
was backwards and up (Y-Z plane, direction of throw). 

The linear accelerations of the CMb for Subject 1 are 
presented in Figure 47. The accelerations of the other 
Bubjects were similar to the one presented in Figure 47 
except about the Z-axis which are presented later in this 
chapter. The acceleration of the CMb became negative about 
all the axes in the double support phase. This indicated 
that in the double support phase the forces were transmitted 
to the hammer which obtained its maximum acceleration at the 


Same time. The acceleration of the CMb in the Z-axis reached 
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maximum and minimum values in different instances for 
different phases and for different subjects. In the first 
turn the maximum occurred after the BDS point while in the 
other turns this maximum occurred before that point. In the 
last double support phase Subject 2 had a continued decrease 
Of acceleration until the release point. This did not occur 
with Subjects 1 and 3, who, after achieving a minimum 
acceleration when the hammer reached its MIN-Z, increased 
the CMb acceleration to a local peak just before the release 
point (Figures 48, 49 and 50). 

Figure 51 represents the moments of inertia of the body 
about the center of mass of the system for Subject 3. About 
the X and Y-axis the moments changed from minimum to maximum 
in every phase. About the Z-axis, Subjects 1 and 2 showed a 
generally lower moment in the first turn than in the other 
urns. In the last turn for Subject 1, in the last two turns 
sor Subject 2 and in the last three turfs for Subject 375 the 
moment became maximum at the beginning of the single support 
phase, then became minimum when the hammer reached the MIN-Z 
point. Increase of the moment of inertia about the Z-axis 
Occurred in the double support during the time interval 
between the MIN-Z point and the BSS point. 

Figures 52, 53 and 54 represent the angular momentum of 
the body about the center of mass of the syste m. There was 
a gradual increase of the angular momentum of the body from 
turn to turn. The global maximum occurred in the second from 


the last double support phase for Subjects 1 and 3, and in 
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X-AXIS(x), Y-AXIS(0), Z-AXIS(z) 
Moments of Inertia of the Body about the CMs, Subject 3 


Figure 51 


tii We ae 


> 
we ; 
fas?) 
° i 
i i 
ve — 
? 
nn 
hep y 
v - 
, = : 
ms é | j 
@ 3 a te ee rem 
4 1 ‘ ad mee . el 
Pad “I -, 1 VE 
z ; PY 
; A's 
Ai. en 
x4 
» 4s a 
y 2 Shy a 
ws 7 tf ees te! 
¥ er, rah are eT dl 7 
5 <4 ae 
~ “a P oad 
d= rad m4 
7 baad ae 
rvy ; 8 eee 
> 4a 
oe - : 
es = —e | vy, 
=f ow am = Lae) ed ere be ee 
~—— = 
7 os . : ad ¥ 
a Ss ‘ 
see an a 
a , & | whe 
Ty, wD ; ‘ 
~ = - ‘ %. 
ia Ww Vv , i. - 
a & ee a — i 
" ne 


- 
{ 


6 


a) 

> 
io 

2 
Ree 


wig 4 


i 
TY. 
iA 
id 

at 


aT 
* 
— > 4 
<= 


ahs 
- 


> 
ee 
ry 


Or 

iG ¥ F — 

SONneF 
1. 


bf 
“ 
ot ¢ 


GR ae ee! 4a 


“~ ~< 
roa ive 
‘ \ 
e 
~ 
a 
en re me TF 
£ 
* 
nd Diag, 
a 
£ 
5 


= 
r) 


a 


eaasseenwive 0 


> 


118 


| 30efqns ‘sWwO 243 yNogGe Apog ©4} JO WN}UEWOW ue NBuY :Z%G aunBLy 


LYOddNS SIONIS 4O ONINNIDSE=Ssg ‘LYOddNS 37gnod Jo ONINNID3S8=Sag 


(99S) JTWIL 
Cot Coo vO°T T8°O 6S°0 


a5 WaT da sdqd 


sui 


O€ 


Sv 


09 


STE 


06 


SOT 


Oct 


Ser 


OST 


W 


0) 


W 


bd 


U 


(Ss /wyw,5y) 


8i1 


4 


-——————— at san 


ass 


RS 


119 


f 30efqns ‘SW 84} jNoGe Apog ay} Jo wn}UsWOW yve_NBuy :€G eunbi 4 


LYOddNS JTISNIS 4O ONINNID|E=Ss¢a “LYOddNS 37gn0a Jo OSNINNI9S398=Sag 
(99S) JWIO 


OC ot LOST 


0z 

OF 
W 
nN 

09 L 
N 
ad 
W 
O 

08 ts 
5 
N 

oot ¥ 

OzT 


(Ss /WyW,b6y) 


4Voe Be 
efi SAPS Ay 


. 
. 
| 


> 
h 
=. @ 
Ee 
i 
i 


a 


¥ a ——— inn a mee fim 
‘2 
= 
rs 
‘ ie € ; 
; Tj 
@s j 
hn ets eR + — - ene: 
ae \ 
_ 
vy } 
i 
vei > 
A 
iy 
ag i 
.. ‘ = as TT <A TH ee . 
- oo } 
‘ = - 
= ’ 
. } = 
' 
' 
4 
i vy 2 ’ ms, 
7 * eg et AE A a ee 2 


" 


ee 


— => = pg ae ——— = 


el ee 


‘et 
i 

. ' 

| 


o 
2 
0 


@ 22 a 
SEGIOW I S MAS 
A*s 
~ Ee 


AS A (oa 


= 
24 

| 

; 


"a 


120 


€ 1y0efaqns ‘SWO 984i }noqGe Apog ey} Jo uN}UEWOW ueLNBuy >7G eunbl 4 


1YOddNS SIONIS 3O ONINNIDSE=SSa ‘LYOddNS 31gnod so ONINNID38=Sdg 


(99S) AWIL 
Get oe oak Sa ue SOL 6L°0 


0z 
W 
OF ‘ 
iL 
N 
a 
09 “ 
0 
W 
08 5 
N 
WV 


00T 


(s/w, Wy By) 
aS Vda 1da sdd oo sdd soa sdqd ssqd sdqd 


aft 


[ 
Lies 2 
: % jl 
| 
| 
| 
| 


| 
| 
| 
| 
i 
. 
4 


a 
4 


r 
~~ 
ag ee 
f 


oe 
— 
rr 


ee —— eee cme + 


rt. 


> 


Le ee i ee A oe a a ie ie ALC IE Ae 


— 


sc | 


oie ~ 


mal 


ae 


= 


Ee 2 > : 


é i) ae 
sd e =. 


- 


121 


the last single support phase for Subject 2. Increase in the 
angular momentum of the body occured during the first two 
double support phases for Subject 1, and during the first 
three double support phases for Subject 3. Subject 2 
presented increases in this parameter during all but the 
last single support phases. Subject 2 also presented greater 
oscillations of the angular momentum during all the throw 
than did the other subjects. During the last double support 
phase, all subjects maintained a relatively constant angular 
momentum until the MIN-Z point. After the MIN-Z the angular 
momentum decreased rapidly until the release point was 
reached. 

The direction angles of the angular momentum vector of 
the body for Subject 1 are presented in Figure 55. These 
angles together with the direction angles of the same vector 
of the hammer (Figure 34) indicate that the planes of 
movement of the hammer and body became parallel in the 
middle of the double support phase just before the hammer 
reached the MIN-Z point. 

The angle of the knee joint of the supporting leg was 
analyzed. During the throw this angle ranged from 75. ts 
m60° for all subjects. In each turn the tangle was greater at 
the BSS point and least at the BDS point (Figures 56, 57 and 
58). All subjects reached the greater angle of the knee in 
the last double support phase. During all but the last turn 
each subject presented a different pattern knee of 


flexion/extention. Subject 1 presented the same flexion and 
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extension in the first and second turn. Subject 2, presented 
a gradual increase in extension after the first turn, while 
Subject 3 presented a gradual decrease of extension in the 
first three turns. Subject 3 started flexing his knee before 
the double support phase had terminated. This could indicate 
that he did not use all the double support time for force 
production. All the subjects started knee extension before 
the BDS point in the last turn. Subjects 1 and 3 reached 
maximum extension before the release point, while Subject 2, 
was still extending his knee at the time of release point. 
The angle of the hip line with the shoulder line ranged 
maemo. tom4/- for Subject. 1, 3° to 50° for, Subject 2eand’6_ 
to 59° for Subject 3 (Figure 59). In each turn the angle was 
greatest at the BDS point. The local minimum angle for each 
turn occurred before the double support phase was terminated 
and the angle started increasing again in the same phase. 
The angles of the trunk with the X, Y and Z-axes 
changed as in Figure 60. Relative to the X and Z-axes the 
angles became maximum during the last double support phase, 
when the athletes applied a final force upon the hammer. In 
the same phase the angle of the trunk with the Y-axis became 
minimum before the release point and started increasing 
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E. Kinematics and Kinetics of the System 

The analysis of the kinematics and kinetics of the 
System was based on the movement of the center of mass of 
the system (CMs) and the movement of the body segments plus 
the hammer about this point. 

Figure 61 represents the spatial coordinates of the 
CMs, for Subject 1. The coordinates of the CMs point for the 
other subjects were similar to the ones presented in Figure 
61. The maximum Z for each turn occurred in the middle of 
the single support phase, while the same maximum for the CMb 
occurred at the end of the double support phase. Figures 62, 
63 and 64 represent the trajectory of the CMs in the 
different planes for Subject 3. Subjects 1 and 2 presented 
trajectories similar to Subject 3. In the horizontal plane 
the CMs moved in an parabolic epicycloid mode. 

Figures 65, 66 and 67 represent the radii of curvature 
of the CMs for all the subjects. The CMs presented the 
greatest translation at the end of the single support and 
beginning of the double support phase. This translation was 
caused mainly by the movement of the CMb, which in this 
instance moved in the direction of the throw. The greatest 
translation among all the turns occurred at the beginning of 
the last double support phase. At this point, Subject 1 had 
greater translation than the other two subjects, achieving 
the highest velocity of his CMb at that point (Figure 44). 

The resultant velocity of the CMs for each subject is 


presented in Figures 68, 69 and 70. The three turns allowed 
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Subject 1 to have smoother velocity of CMs compared with the 
other subjects who used four turns. There was a gradual 
increase in velocity from turn to turn with the global 
maximum in the last turn. The maximum peak for Subject 1 
occurred at the last BDS point. After this point the 
velocity decreased to a local minimum and increased again. 
Subjects 2 and 3, achieved their global maximum close to the 
release point. In the other turns, Subject 1 achieved his 
maximum and minimum velocities at the beginning of the 
double support phase and at the beginning of the single 
Support phase respectively. Subject 2, achieved the maximum 
at the beginning of the single support phase and the minimum 
at the middle point of the same phase. Subject 3 who 
presented the greatest oscillations in velocity, achieved a 
first maximum at the end of the double support and a second 
maximum at the end of the single Support phase. 

The linear velocities about the X, Y and Z-axes for all 
the subjects are presented in Figures 71, 72 and 73. It was 
the velocity in the Y direction which mainly contributed to 
the. resultant. velocity. of..theCMs. There.was.no common 
pattern of movement for the three subjects. 

The acceleration curves of the CMS are presented in 
Figures 74, 75 and 76. The greater acceleration occurred in 
the direction of the throw (Y-axis) for all the Subjects. In 
the last turn, Subject 1 achieved a great acceleration about 
the Z-axis, while in the direction of the Y-axis the 


acceleration was very small or negative. The same did not 
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the X, Y, Z Axes, Subject 1 
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Figure 71 


Lee 


PGE PTET MITE. 
, v2 (Oot ty +a 
. 


4) 
ove 
ty) 4 
’ . cay" 
oh 
Cy 


6 f 
ye ‘ 
ht “ pt 
i ik i ¢ ] 
+ y I 
, reset Ron: Mie 
” se i , 
4 " pn 
Oo eet 
(Vi ba 
: > ae 
“4 7 
i‘) ( 
A 


MING. OF DONETE : 


Je 
pbo*BEGIL 


¥ 


142 


RE LEASE 


BDS 


T IM E (sec) 


BDS=BEGINNING OF DOUBLE SUPPORT, BSS=BEGINNING OF SINGLE SUPPORT 


X-AXIS(x), Y-AXIS(0), Z-AXIS(z) 
Velocity of the CMs in the X, Y, Z Axes, Subject 2 


Figure 72 
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Acceleration of the CMs in the X, Y, Z Axes, Subject 2 


Figure 75 


abr 


~~ 


nN 


QerPepeeruiii 


4 


SDP SBRGIMVINe OF DGABFE 2nb' 


am ry 
Ty 
¢ AM iy 
ees Ny ' 
id Por i 
a | ee S| 
Nh Ee 
ee py wg 
Caletind 


os 


i) 


ot SITY 


ey seer try 
- 


a 
OV entry 


o 
tony 
' ‘ 
Sava WS wh wo Ow - 


eel 


* —— = 
2 Sem oat 

- oe 
: i» 


Se ee ee ne ee 


ey 
a 
° 


RELEASE 


BDS 


Diem cecaces® ° 99 c0 co Fe-on Oo o' 


Fees Qesees 


09 60000000 68 9000 ce Seeteeay, 


«, BBoag rvocesoooet® Bonner" 


T I °M ‘5 ‘(sec) 


146 


BEGINNING OF SINGLE SUPPORT 


X-AXIS(x), Y-AXIS(o), Z-AXIS(z) 


BSS 


BDS=BEGINNING OF DOUBLE SUPPORT, 


Acceleration of the CMs in the X, Y, Z Axes, Subject 3 


Figure 76 
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Moments of Inertia of the System about its Center of Mass, Subject 1 


Figure 77 
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occur with Subjects 2 and 3 who had the same acceleration 
about the Y-axis as about the Z-axis. 

The moments of inertia of the system for Subject 1, are 
presented in Figure 77. The moments of inertia of the system 
for the other subjects are similar to the ones of Subject 1. 
About the X-axis and when the hammer reached its MAX-Z the 
moment of inertia became greater than about the other axes 
in all the turns. After the first turn was completed, the 
moments of inertia about the Y and about the Z-axis changed 
in the same fashion. About the Z-axis the moment of inertia 
became minimum in the middle of the last double support 
phase. This was a reSult of the athletes' effort to apply 
all the forces in the vertical direction which drove the 
system necessarily closer to the vertical axis. 

Figures 78, 79 and 80, represent the instantaneous 
angular velocity of the syste m. There was a gradual 
increase in the angular velocity of the system from turn to 
turn. It became maximum before the last BSS point for 
Subjects 1 and 3 and before the release point for Subject 2. 

Figures 81, 82 and 83, represent the angular momentum 
of the syste m. The angular momentum of the system about the 
Z-axis was Similar for all the subjects. About the Y-axis 
this momentum decreased gradually from turn to’ turn becoming 
minimum in the last BDS point. 

Figure 84 represents the external torques of the system 
for Subject 1. About the Y and Z-axes these torques were 


maximum at the beginning of the throw and were minimimum at 
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re Mw & (sec) 
=BEGINNING OF SINGLE SUPPORT 


BDS=BEGINNING OF DOUBLE SUPPORT, BSS 
Instantaneous Angular Velocity of the System, Subject 2 


Figure 79 
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Instantaneous Angular Velocity of the System, Subject 3 
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Figure 80 
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the end of the throw. The torques remained very small in the 
Single support phases and increased in the double Support 
phases. 

Figures 85, 86 and 87 represent the direction angles of 
the angular momentum vector of the syste m. The angle of 
this vector with the Y and Z-axes started increasing in the 
Single support phase and became minimal at the end of the 
double support for all the subjects. The maximum angles of 
this vector were found to be: 130° with the Y-axis and 41° 
Watbethe Z-axis for. Subject 1, 123° with the Y-axis and’ 35° 
wien the Z-axisyfor Subject 2 and, 125° with the Y-axis sand 
35° with the Z-axis for Subject 3. These angles occurred at 
the last BDS point for Subject 1 and Subject 2, while 
Subject 3 achieved this angle in the middle of the last 
double support phase. 

Figure 88 represents the trajectories of the angular 
momentum vector for different subjects. These figures 
indicate that while the system rotated about an axis which 
was defined by the angular momentum vector, the angular 
momentum vector was also rotating about another axis. This 
axis passed through the CMs point and intersected the ground 
in the center of curvature of the orbit described by the 
ground intersection of the angular momentum vector. The 
Second axis was the axis of precession of the syste m. 

The total kinetic energy of the system for Subject 1 is 
presented in Figure 89. The energy was mainly rotational 


rather than translational (Figure 90) and inreased gradually 
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f the System, Subject 1 
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in the same way as the angular momentum of the Syste m. 


F. Discussion of Hammer Throwing 

The velocity of the hammer increased gradually from 
turn to turn, with a local increase in the double support 
phase and a local decrease in the single support phase. The 
global maximum velocity occurred 10 to 20 milliseconds 
before the release point. Theoretically, the hammer should 
have reached maximum velocity at the instant of release. The 
facts that*’atedid not could! be a result of data smoothing, of 
misdetection of the release frame or that the analyzed 
subjects did not utilize the maximum velocity of the hammer. 
To ensure that the timing of the maximum velocity was not 
affected by the smoothing of the data, the velocity was 
recalculated using raw data. This test confirmed the time of 
occurrence of the maximum velocity to be as presented above. 
By filming at 100f/s, at the release point the hammer moved 
= 30°cm from one frame to the other. If there was a 
misdetection of the release frame then the hammer had been 
released in an interval of =~ 60 cm. However, in such a case 
the height of release would have been = 80 cm which seems 
unreasonable to believe. Also, assuming negligible air 
resistance during the flight of the hammer and by using 
equation (1), the calculated heights, angles and velocities 
of release could produce throws over 90 m. It is therefore 
reasonable to accept that the athletes did not utilize the 


maximum velocity of the hammer at the release point. 
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The height of release was lower than the standing 
height of the athletes' shoulders and it depended on the 
acceleration of the hammer in the vertical and in the 
throwing directions rather than on the height of the 
athletes. 

The angle of release, although none of the athletes 
achieved the calculated optimal, was a critical factor for 
the final throwing distance. To improve the angle an athlete 
must learn how to break his horizontal motion and also how 
to achieve an angle of the plane of motion of the hammer 
which is closer to the optimal angle in the last turn. 
However, the length of the hammer might prevent this angle 
being achieved during the last turn since the ground would 
interfere. 

A faster single support than double support phase for 
all the turns was a common characteristic of the three 
subjects. This was achieved by using the last part of the 
double support phase for initiating a fast turn. During the 
first half of the double support phases the CMb had achieved 
the greatest velocity in the direction of the throw. The 
torque of the CMb about the supporting foot created by the 
new position of the CMb with respect to the supporting foot 
obliged athletes to make a fast landing of the free foot. 
The free foot landed behind the supporting foot in the (Y, 
-X) section of the circle and in an angle of 70° measured 


from the Y-axis in the direction of rotation. 


ad 1h ae es 
RAEN bs1% ; 
Vel “ 


oe 7 
: 


Ps ba 


piibneds <i 0 rawol cas senator 3 
add Ao bebneqeh $i Snes suebiueds 395 jotta of, i 
oa? at bas Isottvev odd: ak setamet aa 20 8083 
odd do Adgist hlie'te eal yedtes oooh 


s i 
y ' 


B 
he 


sesoidts add Ie, anon souorts te sessiex 190 efor 


i ee 


160% yosael tentatrve reer _ Lamidqo” bataltiate>:, edt? 6: 


hs 


sisifis 6 signe, ems avorgmi oF ‘anneal Rohe EH 
won oals bre nektom ininos iyod eid aaerd ot worn 
remmed aff fel m@itom io sasia, ond 35 elpas ne iste ) 
vue Beek get ai sions Lemiaqdiens od 320k: pe 
elyonw site geaweag ‘tap bn Tema adhe dipnel oft 
biiueie Rang ast agente wuts “deat edz. enti bevatits 7 
19? semig S40GRve Sldvod vada s2:0qque slonie 1 , : 
saxda a03 29 chagiies2a1ehy nome? A 26% ue 
ait togase 2460 ot prioy a (Bavaidas aaw aad? 
a3 Dal jwiaat tes? & gpnivetsing as, szarig rn 
bsveitss bei: fs, agedcry 220¢que eLduob of to + 
ont iil ede, Id noiiseri®. sda iat x3 toolayigesda 
ald yd Beiseys tn02. pslozcqquel sd guods dM3 eat a0 
ck sdcadaneiill bist bp Josgees difiiw dO edt 20 ‘nobdiae 
Seek aay Re eh pete fast @ etsy 02 estold: 
a add ni d6o2 eatstoqqua silg(Badded bebaat toot oath 
bewizene °OT 36 sivas oa Ab BAgystoais att 1oinotoan (| 
Mofteiox 26 eine tens otae-t, adistont 


i 4 f ; hs \s > i 


is) 


The acceleration of the hammer increased gradually from 
turn to turn with maximum local peaks in the double support 
phases. In the single support phases, athletes accelerated 
the hammer by rotating very fast and by reducing the moment 
of inertia of their body about the vertical axis. It was the 
duration of acceleration in the double support phase rather 
than the maximum acceleration that improved the final 
throwing distance. 

The shoulders were continuously driven by the hips, 
while the hammer was continuously driven by the arms. This 
was accomplished by a very fast legs/hips rotation in the 
Single support phase which allowed the athletes to have the 
greatest possible hip-shoulder angle at the beginning of the 
double support. It was only at the changing points from one 
Support to the other that the hammer assisted the rotation 
of the system slightly. 

The trunk was utilized by the athletes to transfer the 
angular momentum from legs/hips to the hammer. It received 
angular momentum from the legs/hips in the double support 
phase and transfered it to the hammer in the Single support 
phase. Only in the last double support phase, after the 
force production from the lower limbs had been terminated, 
did athletes use the trunk as a force production system to 
improve the final velocity of the hammer. After this point, 


the shoulder musculature participated in increasing this 


parameter. 
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The knee joints achieved maximum extension at the 
release point, while in the other turns extension was less 
than maximu m. This knee flexion/extension pattern helped 
athletes to rotate faster in the single support phases. 

The moments of inertia of the body about the X and Y 
axes reached a local maximum and a local minimum in every 
Phasevotvthe throw. About the vertical axis and’at the end 
of each double support phase the moments became maximum and 
decreased from there to the end of the double support 
phases. This helped athletes accelerate the hammer in the 
Single support phase, even without utilizing any force 
initiated by the legs. 

The hammer and the CMb rotated in parallel planes only 
at the point before the hammer reached its lowest point in 
the double support phase. 

The system:sathlete + hammer translated at the end of 
the single support and beginning of the double support 
phase® This®translation was greater for the*three turn 
throw, because the athlete had more freedom to move. 
Translation in the direction of the throw was least in the 
last turn for the three turn throw. The fourth turn did not 
allow athletes to break the translation of the system, 
resulting in a lower angle of release and probably a lower 
release velocity of the hammer. 

The angular momentum of the system increased in the 


double support phase and decreased or remained relatively 


constant in the single support phase. 
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The direction of the angular momentum vector of the 
body came closer to the vertical at the middle point of the 
Single support phase. At the same time the direction of the 
angular momentum vector of the hammer was at its greatest 
angular displacement from the vertical axis. 

The torques for the rotation were mainly applied to the 
System during the double support phase. 

The maximum angular velocity of the system for each 
turn occurred at the end of the double support phases. The 
time of occurrence of the global maximum angular velocity is 
important because it allows athletes to transfer the angular 
momentum of the system to linear momentum of the hammer. The 
result of bad timing was reflected in the radii of curvature 
of the hammer; Subjects 1 and 3 achieved a greater length of 
badricthat «did Subject )2. 

A fourth turn served basically to smooth initiation of 
the throw rather than for increasing the dynamic parameters 
Diathe syste ms This turn could be detrimental in final 
distance since at the release point, athletes had the 
tendency not to break the horizontal movement. In the last 
turn, linear momentum of the system could not be transferred 
to the hammer. The cost could be a reduced angle of release 
of the hammer and possibly a reduced release velocity. A 
common mistake for Subjects 2 and 3 was that their CMs 
reached maximum velocity just before the release point, 


while Subject 1 achieved this maximum at the beginning of 


the last double support. 
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A heavy athlete will have some advantage over a lighter 
athlete although this does not seem to be a critical factor 
Since a lighter athlete can balance the lack of great mass 
by rotating faster. The height of the athlete does not seem 
to be a critical factor either. Although a tall athlete 
could theoretically achieve a greater radius of rotation of 
the hammer and a better angle of release, the present 
analysis has demostrated that this did not happen. An ideal 
hammer thrower would have very strong and fast legs, very 
good coordination, and a strong lower back. "A strong kid 
who iS a very good dancer". 

The system:athlete + hammer can be characterized as a 
heavy assymetrical top which describes precessional motion. 
Since the angle of precession of the system varied, another 
characteristic of the top is nutation. The axis of 
precession of the system can be used in the coaches language 


acethe (axis eofMrotation. 
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V. SUMMARY AND CONCLUSIONS 
The purpose of the present Study was to develop a technique 
for analysis of the hammer throwing event and to analyze the 
technique of world caliber hammer throwers. To develop the 
technique for analysis, the DLT 3-D cinematographic method 
was tested. These tests included different mathematical 
models for image refinement, as well as the use of different 
calibration trees for the derivation of the calibration 
coefficients. It was found that a tree, geometrical in 
shape, together with mathematical model IV for image 
refinement (Karara and Abdel Aziz 1974), gave the most 
accurate results. The RMS error of meaSured versus simulated 
coordinates of the control points was: RMS(X) = 0.24cm, 
RMS(Y) = 0.19cm and RMS(Z) = 0.26cm. Six points were chosen 
to be treated as unknowns in the outside area of the 
calibration tree. The error of these points was found to be: 
RMS(X) = 0.27cm, RMS(Y) = 0.60cm and RMS(Z) = 0.72cm. The 
latter error indicates that a smaller calibration tree can 
be used to calibrate larger areas. 

The calibration tree which gave the best results was 
used for the calibration of two cameras in the data 
collection of the hammer throwing study. The data were 
collected during the 1982 European Championship and the best 
throw of each of the three medalists was analysed. Two phase 
locked cameras were used for data collection. A Lagrangian 
polynomial was used for the time matching of frames, based 


on an external event which was a light of 10 Hz generated by 
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an external light generator and recorded on the sides of the 
films. The basic DLT equations together with model IV for 
image refinement (Karara and Abdel Aziz, 1974) were used for 
the transformation of the digitized to real life coordinates 
of selected body and object points. Digital filters were 
used to reduce the noise inherent to the data. The derived 
formulas for analysis of the data were based on rigid body 
dynamics and can be used for analysis of any rotational 
movement. 

The results of the analysis of hammer throwing can be 
summerized as follows. 

1. A faster single support phase than double support phase 
waS a common characteristic for the three analyzed 
Subjects. The torques of the CMs with respect to the 
Supporting foot were the main reasons for the fast 
Single support phase. 

2. The height of release was lower than the height of the 
shoulders and it depended on the acceleration of the 
hammer in the vertical direction rather than on the 
height of the athletes. 

3. The angle of release was a Gr pthtica hifactomaior rrhe thina i 
throwing distance. However, the length of the hammer 
might prevent this angle being achieved during the last 
turn since the ground would interfere. 

4. The angle of the plane of rotation of the hammer started 
from an acute angle with the horizontal plane and 


gradually changed from turn to turn to the release 
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angle. 

The length of the radius of curvature of the hammer 
decreased from turn to turn until the instance when the 
hammer reached the minimum point in the Z-axis during 
the last turn. After that point the radius started 
increasing to become infinitive at the release point. 
The hammer was continuously driven by the athletes and 
it was only at the changing points from one support to 
the other that the hammer slightly assisted the rotation 
of the system. 

The velocity of the hammer increased gradually from turn 
to turn, with a local increase in the double support 
phase and a local decrease in the single support phase. 
The maximum velocity occurred 10 to 20 milliseconds 
before the release point. 

The acceleration of the hammer increased gradually from 
turn *to turn and in the double support phase. In the 
Single support phase, the athletes improved the 
acceleration of the hammer by reducing the moment of 
inertia of their body about the vertical axis. 

The hammer and the center of mass of the body rotated in 
parallel planes only at the point before the hammer 
reached its lowest point in the double support phase. 
The maximum vertical displacement of the CMb in each 
turn was less than 0.12 cm. 

The radius of curvature of the CMb had noisier 


oscillations in four turn throws than in a three turn 
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throw. 

The linear velocity of the CMb increased during the 
Single support phase in all the turns. 

The moments of inertia of the body about the vertical 
axis increased during the second half of each double 
Support phase. During the single support phase this 
Parameter decreased continuously until approximately the 
middle point of the double support phase. 

The angular momentum of the body increased from turn to 
turn reaching the maximum during the second from the 
last double support phase. 

The knee joints achieved maximum extension at release 
point, while in the other turns extension was less. Only 
in the last double support phase, after the forces from 
the lower limbs had been transferred to the hammer, did 
athletes actively use the trunk as a force production 
system. 

In the horizontal plane the CMs moved in a parabolic 
epicycloid mode. 

The system: athletethammer translated at the end of the 
single support and beginning of double support phase. 
This translation was greater for the three turn throw, 
because the athlete had more freedom to move. 
Translation in the direction of the throw was least in 
the last turn for the three turn throw. A fourth turn 
did not allow athletes to break the translation of the 


system, resulting in a lower angle of release and 
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probably a lower release velocity of the hammer. 

There was no common pattern of linear velocity of the 
System among the analysed subjects. 

The vertical acceleration of the system during the last 
turn and the minimum horizontal acceleration in the 
direction of the throw in the same turn seem to be 
Critical for a successful throw. 

The maximum angular velocity of the system for each turn 
occurred at the end of the double support phase. 

The angular momentum of the system increased at the end 
of the double support phases and decreased or remained 
relatively constant in the single support phases. It was 
the angular momentum of the body that directed the 
motion. The direction of the angular momentum vector of 
the body came close to the vertical at the middle point 
of the single support phase. At this time, the body 
rotated mainly about an almost vertical axis. During the 
Same time the direction of the angular momentum vector 
of the hammer was at its furthest angular displacement 
from the vertical. 

The torques for the rotation were mainly applied upon 
the system in the double support phase. 

The kinetic energy of the system increased gradually 
from turn to turn and it was mainly caused by the 


rotation rathet than the translation of the system. 
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On the basis of the results of the hammer throwing analysis 


the following conclusions seem justified. 


i 


An athlete should try to reduce the time of the Single 
Support as much as possible. This can be achieved by a) 
tilting backwards during the first half of the double 
Support; b) by using the second half of the double 
Support phase for a fast rotation; c) by rotating the 
Preeelegcclosecto othe supporting! legrcd)isbyeplacing the 
free foot behind and on the right of the supporting foot 
and at a angle about 250° to the direction of the 
throwing axis (measuring the angle from the back of the 
circle and in the direction of the rotation). 

A fourth turn is basically used for a smooth initiation 
of the throw rather than for increasing the dynamic 
parameters of the system. This turn might cost in final 
distance since at the release point, athletes have the 
tendency not to break the horizontal movement. The cost 
could be a reduced angle of release and possibly a 
reduced release velocity of the hammer. 

It is the duration of acceleration in the double support 
phase rather than the maximum acceleration that improves 
the final throwing distance. 

A heavy athlete will have some advantage over a lighter 
athlete, although this does not seem to besatcrtihical 
factor, since a lighter athlete can offset the lack of 


great mass by rotating faster. The height of the athlete 


does not seem to be a critical factor either. 
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5. The length of the hammer might prevent athletes from 
achieving an optimal angle of release since the ground 
would interfere. Therefore, this length might have to be 
reduced in order for athletes to perform mechanically 
efficient throws. 

6. The system: athletet+hammer can be characterized as a 
heavy assymetrical spinning top which describes 
precessional motion. The angle of precession of the 
System varied, therefore the system had another 
characteristic: nutation. 

7. An ideal hammer thrower could be an athlete who has very 
strong and fast legs, who has very good coordination, 
who has a strong lower back muscle system and, who is 
not necessarely very tall. "A strong kid who is a very 


good dancer". 


Recommendations 


The following studies are recommendations on the basis 
of the present findings: 

1. The DLT method should be tested on a statistical base to 
estimate the expected error in areas which are not 
covered by the calibration tree. 

2. New methods should be found for the calculation of the 
anthropometric data of the body segments, based on 
individual characterisics of the analysed subjects. 

3. A statistical study is recommended for comparison of the 


kinematic and kinetic parameters of different level 
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hammer throwers. 

Direct force measurements should be used for the study of 
the forces involved in the different phases of the throw 
and their direction. 

A study of hammer throwing in the last turn of the 
movement with a filming rate of at least 300 f/s is 
recommended to examine whether the highest velocity of 
the hammer does occur before the release point as it was 
found in the present study. If the results of such a 
Study agree with those in the present study, then 
biomechanists, coaches and athletes should seek new 
throwing technique which will allow athletes to utilize 


this maximum velocity at the release point. 
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Description of the DLT Technique. 

This description is based on the technique invented by 
Abdel Aziz Karara in 1971 and presented by other 
investigators (Marzan et al, 1975; Walton, 1981: Miller et 
al, 1980). The technique is based on the basic theoretical 
concept of photogrammetry, that is: the photograph, being a 
perfect plane, is a central projection of the object space. 

Figure A.1, iS a geometric representation of the 
transformation of the 3-D object space onto the 2-D 
secondary image. The following manipulations are based on 


vector analysis. 


Let: 
Pterxkay Toynpwsizn. (AO) 
Q = an, + Dnm + CNn (An) 
R = XoNx + YoNy + ZoNz (A.3) 


Nhere-icx,y,2,a,b;C)Xor¥orZo, are Scalars 


andenje(i=xyy,z,k/m,n) unit vectors. 


If s = the magnified principal distance then: 


S =a-sn, (A.4) 

Let: 
Pp = un, ae WAR: (A.5) 
(A.6) 
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where: u,vV,Uo,Vo are scalars 
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From figure 1 we have: 


By Pes Sg (ATI) 


meomevA.5)¥ (A .'6) ,CAs17) 


Ges GU ag) net ye (A.8) 


Further observe that: 


LN Jghce audake (As. 9) 


From (A.4);,(A.8),(A.9) 


es (UPS Aig.) Net SOF Sey co ee Sn, (A210) 


Assume that n, and n, are orthogonal; then equation (A.10) 


can be written as follows; 


Tee GUa\a sy POLY evo) Ngee Onl CA amid 


sing equations: (Az 11)~and (A.2),; and’ thescolinearity 


conditions, we obtain: 
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as 

a oe (A.14) 
cS 
bs 

Niep Vege = (A.15) 
Cc 


u and v are the scalar components of the vector p which 
locates the image point I, with respect to the image 
reference frame. In digitizing, we have digitizer units not 
necessarily the same as those used to establish the object 
reference frame. 

Pntrcoduce scale factors C,, €z for n, and n, 
directions, respectively, and let; 

U = Cu (A. 16) 


V = C2v CARE) 


Brom equations (A.14), (A.15), (A.16), (A.17), 


as 
U-= 3G; Sue - —) (A. 18) 
eS 
bs 
Vesechiveents } (A.19) 
G 
Equations (A.18) and (A.19) express the digitizer 
coordinates U,V in terms of a,b,c. 
From Figure 1, 
(A.20) 
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Promeequationsi,(A.4)%..¢A42),.¢(a.20): 
O-= tax-teblay t+ tyja polny +2 (z.— zodnz Chea) 
Let 61,,, 912, 6,3, be the direction cosines of ny with 
respect to the object space frame (ny, Mm, Mn) } O21, O22, 
623 be the direction cosines of ny with respect to the same 
Erame; and!0@3; @32,°033,-be°the direction cosines of ngaw ith 


respect to the same frame. Then: 


Nx = O14Ne + Gy2Mm + )O43Np (As22) 

Ny = G241ne + G22Nm + G23Nn GAs 2 3)) 

Nz, = O31N, + O32Nm + O35Nn (A.24) 
mae 0; 5°(151,2,3,(j5?,2,3)eform/a 3-dimensional rotation 


Matrix. 


From equations (A.21), (A.22), (A.23) and (A.24), we obtain: 
Cee Blk x= Xo) 04 PoP by = VolOpe hf A] = 26 )0 sane 
+ C(xbhs *o}@sdet Cy - yo)@22 + (Zz —e2ig)) 052m 


Be Be 7 x@) 045 & (y va Yo)823 5 At > ZovGss NG (A.25) 


Equate the coefficients of equations (A.2) and (A; 25) 


ae sk (x Ko) Cai + (y eS Yo)@a1 til oars Zo) 431 (A226) 
b =Bdxarixok0i2.e9lyan Yo)922 + (Z =)Z0)O32 (As27) 
ensfixmatxo)@,4,0tT lyes ye lessor (2 - Zo)633 (Ae 23) 


From equations (A.14), (Agt5den (Ae26)eesAa2eheand. (A.28): 


§ et 

7 
Ace a 

: a 


uaa 4 Rata 


x8! ; + b a 
. See bb al os ve 
2 (08; cA) ose) 


(1.0) antes + a) # net a ea 


d3iw\.q to senkeeo moitosxtb ads dey 


seo ce? 34 wall val) omard rv rs th 


suse edd o3, tosqeet dstiw yn. 20 eenteoo eee Te 
1} etre? "00 ee eat aS x 


dtiw, »a 20 eani@o’ otsaerit es od en ae 
\ 4 a # ane PF . ; 


nett,» i i ' ig 
2 : } ey Bigs 7 
(ss A) ’ \ afle c < ae ay a : 
a eT | t 
(ES.A) nfes® = att 4 + said = et. 
(85.4) wes 94 tear. 
noitedos tgubidientt~ és w30% 30x, “t 8, 
ead ro +i 
J + . ee 
Ab pee. Fires 7 
inisido ew — brs Ros cotta ae : 
: ~‘T ies we ba 


hia hei +) on the 2 ae 
ailveties > 3) + satin = 
(25.4)  ppteeWGet) > A) + cstlak = 


ASS a 


(aga) Baa US.A) sie er 


(iS) 
(TS A) 
(es.8) 


where: 


with: 


L,x Ug L2y + L32Z de ie 
So ee REE A aeRO! Ok Teak 
Lexet) baevitoLaiztse 1 


Beko Ley t bye Hike 


Lexetsbaoy taL4qzet+ 1 


Ly wetG, ¢ugOqsb-ts0%94/a 

EZ =h@,4ug6ederss0e,)/d 

ReteeGy (us@s$uo $5407 

Bes=tneCa [itued. 4. anseores ite (goss 
tobuge3 sncisés,dzoiva 

Bs e=nhC2 (vo samned 15)/d 

Bad=fCoivebe; vaises.07d 

LyverCetvediser séssiy/d 

Lyos cei bivieees. Helaudxe tolveeas 
ti(veesjaz s@oedzodyad 

Be 2céas7d 

Heoo= 6439/0 

Bye@n=66537d 

Gd = -(xX0613 + Yo923 + 20433) 


188 


(A.29) 


(A.30) 


(At8a9 
(A.32) 


(As33) 


~ SO21)Yo 


(Ai34) 
(A.35) 
(A.36) 


(A237) 


= SO022)Yo 


(A.38) 
(A.39) 
(A.40) 


(ASA) 


(x42) 


Equations (A.29) and (A.30) are the basic Direct Linear 


Transformations (DLT) equations. 


L, (i=1,2,.....11) are the transformation coefficients 


expressed in terms of seventeen separate elements. 


(Xo,Yo,Z0) define the location of the camera with respect to 


(OS.A? 


(O&A) 


-=% rama coe 1 


(tea) | ; dante "BC th eB e® Bie: 
(SE. A) | | | EM ae wikabep hed = 


(£€.A) : | | BY ¢ he: pia 
ov (1 e828 = ee + ght) 8p “otb)lia> = 
(BELA) | \ Lab Be, th eee: * oo | 
(ab lay) to A) Seated ~ Fibov) wi} 
(ae A) DAL B\ (e_82 y Mod lal 
(tea) 's uN ay 6\ tg etanmiaa hd dca as 
ottsete > eto . ot (sre ag Pony ei ie ai Bop 
(HEA) er, BN os aah ATasBovdow\ 1h i 
(@E.A) | 20 al ee fio BNeb b+ wy , 
(08 A) “St tilde Maa ‘i PBN cab moval 
(1e.A) re. 7 me B\ec8 = al l rm 
r ti = ie ote i hai 5 
($0.4) | teetos * onto * cyton) te 


tss0il toe 72?C siesd er : * 
nga 
2jnsici3iecs Aotiemvolenend 


62 toege02 dain ovemss ent? 


189 


the object reference frame. 

(Uo,Vo) define the location of the Principal point. in the 
secondary image. 

(C,,C2) define the scale factors. 

S defines the magnified principal distance. 

9;; are the directional cosines defining the orientation of 
the camera, with respect to the object space. If the vectors 
n, and n; of the secondary image are not orthogonal, then 
the angle between the two vectors is the eighteenth element 
which express the eleven transformation coefficients. A 
minimum of six control points is required in order to solve 
a system of 2(K) [K = number of control points] linear 
equations and find the values of the eleven coefficients. 
This is an overdetermined system and to solve it a least 
Squares method can be used. However, to reconstruct the 
three dimensional spatial coordinates, at least two 
observations from different points are required. Here two or 
more systems of 2(K) linear equations each have to be solved 
for the eleven coefficients of each camera. The DLT 
equations then are used to solve for the unknown (x7 Vy so 
spatial coordinates of points in the space. In this case a 
system of 2(J) [J = number of observations] linear equations 


has to be solved for three unknowns. The same least squares 


method can be used. 
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Figure A.1 Geometrical representation of the transformation 
of the 3-D object space onto the projected secondary image. 


A = Origin in the object space; 
B = Origin in the image plane; 
N = Projection center; 


Mp(uo,Vo) = Principal point; 

Oix,y,z)"= A point in the object space; 

I(u,v) = Photographic image of point O; 

kr7Om,;N,» = Right handed orthogonal triad; 

ny,N, = Right handed orthogonal triad of unit vectors; 
nN, = Image plane orthogonal unit vectors; 

= Position vectors defining relative positions of 


Position vectors defining relative positions of 
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Solution of the Linear System 

Let N be the number of the CP used then the system of 
equations takes the following matrix form: (SOR WEED are 
the real life coordinates of the CP (i=1 to N number of 


C.P); U;, V; are the digitized coordinates of the same 


points.) 
meen) OURO G08 40 =XGUy rE P1472 20), U, 
Pome fee?) OVEQe*OL sO —X9Us p-YSU 7 2250s a U2 
mei Wh kU OQ) 60 Oe -X2U5 -¥5U5 -Z25U3) 15 U3 
SESE Teter ate eve! “oe ia wi lie Gs 60) s/s 6 oe jelele ste Me eters | tre 
REGUS EMC C! sist 6 6 cies e.s ave celeleve's) eleietele: | eters sami iiae ats 
= (A.43) 
emer cot) OO. 0 SAnUn Tipe nee one na tis UF 
Ome O70 Fy OYA pel sR yVa SVG Vee S257 eras V; 
Dae mrOrwO ek V5 ° Zo 1 “XSoVo5 7 V2V oe -25NV 2) [Ge V2 
Ome OarO: “OCX PVs Zs) HXoV3 KY sVGysSsVaibs V3 
a stette es se ae 6 ° gies ate tte 6 ° Saba een we aS 
DeeuO mer ae Kee Toe tT Ry Minn oi hy nea nel oan U, 
[Ad (bel stece 


The solution of the above system is given by the following 


matrix equation: 
B= (A'*A) @ °C) (A.44) 


Z rs | 
where: A' is the transpose of the Jacobian matrix [A] and (A °A) is 


‘ 
the inverse of the (A -A) matrix. 


The same equations are used for the second camera. After the 
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Mee PEO etood 1) sandl thes! +o(tbe1lté 11) PoEt thesfarstmandt the 
second camera respectively have been calculated, the 
reconstruction of the 3-D real life coordinates of the CP is 


achieved by the following matrix equation. 


CUES ¢ 7 Eve) (UE4 o> Lite) (UL4 ¥> i) Cla a U) 
x 
OVES = aah) (VLEe= ney LVL 4S ee Ga ae Vv) 
vie (A.45) 
(UL = L,) CU +67 As) (OLY ,- a) (Li ra U) 
Z 
(VE. 5 ev (VEdo> ig} (VL44- bess) Cie ee V) 


[E] [F] [G] 
the solution then is given by the following matrix equation: 
E-'°G = F (A.46) 


where: E~' is the inverse of the E matrix. 
In the equations of the system (A.45) the L,; (i=1 to 11) and 
U, V are the coefficients and the digitized coordinates of 


the CP from the first camera and 1; (i=1 to 11) U, V the 


Same for the second camera. 


Image Refinement 


The coefficients of the equations (A.29) and (A.30) 
count for the transformation of the 3-D object space onto 
the 2-D secondary image and vise versa. However, in fiim 
analysis the above coefficients do not count for the 


distortions of the image due to the lens of the cameras and 
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ENED ojectoreas iwellcasstoithe film itself) (Werassumed that 

the film is a perfect plane which is not always the case.) 
For the refinement of the image, Abdel Aziz and Karara 

(1974), have presented 6 mathematical models to be used 


together with the DLT equations. These models are: 


MODEL I: Linear polynomial in U and V: 
OUs= a ;etVas0 * av CA age) 
AV =saage teavUabiacv (A.48) 
In this model, only the linear components of lens distortion 
and film deformation are taken into consideration. 
Incorporating Model I into the basic DLT equations does not 
change their form. In other words, by uSing Model I in the 


DLT solution the number of unknowns remains 11. 


MODEL II: One more unknown k, 18S added to Model I to account 


for symmetrical lens distortion: 
AU = Qtr a2U ay a3V = Rkank (A.49) 
AV*= ay t+ @eU +7 acV FYE ee (A.50) 


In this case, the DLT solution involves 12 unknowns. 


MODEL III: An odd radial polynomial of the seventh degree is 
added to Model I to account for symmetrical lens distortion: 
AU = eyRthasU + GV 't+Xtkyr? ot kkye Steceey) A651) 


AV = a,0+(a@2Uet2aéV * Vbker? + -kor* Soker) (A 62) 


The DLT solution in this instance involves 14 unknowns. 
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MODEL IV: Model III is combined with a model for 
asymmetrical lens distortion: 
BURKS a4 2 a2U" + asV +O (kart ek yee ee, 
+a PON MEGS) + 2P> xy (A 53) 
AVLSGay Fla, USt. AGVEFLV(KEE2OF KOE ‘DSIRE red 
tEPU tre $8272) “hePi my (A.54) 


The DLT solution here involves 16 unknowns. 


MODEL V: The same as Model IV, except that the radial 
polynomial accounting for lens distortion is a full 
polynomial of the seventh degree: 
AUTSsilagrt® Gourtiasyv 
ere ( Koee +c keres+or. eyerery ares Kerns) 
+ oPetr out (2% 2) ky (A, 55) 
AVGela,CteacU foaEgVv 
£ ¥Ck, peOvrekMPeecak }escepe, MIVGGL ks rss 
+ P(r? + 2¥?) + 2P,XY (A.56) 


Here the DLT solution involves 18 unknowns. 


MODEL VI: Same as Model V, except the polynomial in U and Vv 
is of the second degree, accounting for the nonlinear 
component of film deformation: 
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t Ber topr2Va)cne2Pexy (A.58) 


In this model, the DLT solution involves 22 unknowns. 


In all the above models 
U and V are the digitized coordinates of the point under 
consideration. If Us & Vs are the image coordinates of the 
Be tnt of symmetry, then 
BherUaseUssandsY =nVr=.Vse 
r = the length of the vector from the point of symmetry to 
the image point under consideration. 
AU and AV are image refinement components. 
meiiaiy...,10) are tcoeffichetitsaof film deformation. 
k; (i=1,...,5) are coefficients of symmetrical lens 
aistortion. 
P,, P. are coefficients of asymmetrical lens distortion. 
When one of the above models (except model I) has to be 
used with the DLT equations, then an over-determined systems 
of 2(K) non-linear equations have to be solved for the 
appropriate number of unknowns for each model. Introducing 


one of the models for image refinement, equations (A.29) and 


(A.30) take the following form: 
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where: X,Y,Z are the coordinates of the point in space. 
U,V are the digitized coordinates of the point. 


Beate 5%, cit) are coefficients named transformation 
coefficients. 


AU,AV are image refinement components of lens 
distortion and film deformation. 


Using one of the mathematical models for image 
refinement (exept model I), the number of the coefficients 
increases and an increased number of CP is required. Table 
A.1 contains the number of coefficients and the minimum 


number of CP required for each model. 


TABLE A. 1 
NUMBER OF UNKNOWNS AND MINIMAL CONTROL POINTS 


FOR EACH MATHEMATICAL MODEL 


Model Unknowns Control points 
No Ref 11 6 
Model II 12 6 
Model III 14 7 
Model IV 16 8 
Model V 18 9 
Model VI ge ‘a 


For the optimal number of the object-space control 


points, a mathematical model is given by Karara and Abdel 


Aziz (1974): 
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where : S, is the standard deviation of the standard 
deviations of the object-space coordinates 
S is the standard deviation of the object-space 
coordinates. 


n is the number of observations (in case of two 
cameras, 1t 1S twice the number of control points ). 
u 1s the number of unknowns. 


The same investigators have indicated that beyond 20 to 
25 CP the improvements in the results are probably not worth 
the additional effort required to position the extra control 


points. 


Solution of the Non-Linear System 
Assume that model II is to be used. The equations 


(A.59) and (A.60) can be expressed in the following form: 


F, Pigghe ts 5 ths Z + thy 

nt oka kr (Lk + Ly OW te banZ Fie = 20 (A.62) 
F> Site me erase Vart Ib 21% elegy 

PCy + ki Yr?) (Gok + Laghgt tare * ee (A.63) 


Note: Model II doesn't include a;, a2, @3, Gy, as since 


these are components of the linear lens distortion and film 
deformation which are included in the basic 11 coefficients. 
By expressing (A.62) and (An63) in “truncated Tayloreseries, 


we obtain: 
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Baby tAL;, L2tAL2, LatAL;, BetAL,, LetAbs, LetAbe, bz tAL7,; 
Sproles, Lethbe, LyotALyo, LFetAL,;,, ky+Ak,) 
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ell 4) tess 0's. (A.64) 
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Let J be the above matrix, then if G is the 
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where J° is the transpose of the Jacobian matrix. 

The elements of the C matrix and in the first iteration are 
calculated from equations (A.62) and (A.63), by inserting 
the initial values of the 11 basic coefficients as have been 
calculated from the system of equations (A.44). Initial 
value for the k, element is used the zero value. The 
computer program is set to attenuate the iterations when the 
sum of the absolute values of the residuals is smaller or 
Peedgivro the criteria value of 10°*°;: (21Ri = 10>") stor 


when the system stops converging. 


Simulation of the Spatial Coordinates 

After the choice of the appropriate model for the image 
refinement a number of 2U system of equations is solved for 
the reconstruction of the X, Y, Z coordinates of the points 


under considaration. 
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equations: 


Libg + Lelio + L3L;, 
awa. se (A.80) 
Eig ft denoe) Tela 
LsbLg + Lelio + L7L;, 
NE a aca ceri (A.81) 
Lig tt Lie, abot kei 
The same way for U, and V, of the second camera. However, 
the above equations are not valid if the CP are not 
Symmetrically distributed in space which is usually the 
case. Therefore, one must enter the coordinates of the 
center of symmetry manually. 
After the calibration of the system and the digitizing 
of the coordinates of the points under observation is 


finished, the system of equations (A.79) is used to simulate 


the spatial coordinates of each point. 
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DLT Programs 


The following programs and subroutines were written in 


HPL languege for the H.P.9825B desk top microcomputer. 


Se No a ee eee eee 
DATA ENTERING PROGRAM 

Tape D.Lur sik. 

track #0 


file #0 


DESCRIPTION . 
This program was designed to insert all the required 
information for the DLT calibration. 


User enters: 

be elrhe = of “all CP; 
2. The names of all CP; 

3. The X, Y, Z measured coordinates of all CP: 

4, User and project I.D. information. 

User digitizes: 

1. The U, V coordinates of all CP from both cameras; 

2. The coordinates of the edges of each film for the Us and 


Vs. 


The computer prints all information (except users and 
projects I.D. If you wish this information in print, it will 
be recorded in trk 1, file 0, string A$[10, 60]). 

After each entry the program is set to ask whether or not 
you have entered a parameter with error. If yes, then you 
should type in Yes or Y or y, and enter the correct 
Darameter. This correction feature is not in use for names 


Of the CP or 1.D. information. 
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% WHEKKKKRKEKERKKEKKEKEEKKD OT TOI CK A TREK O 


eo -D CALIBRATION (DATA ENTERING)"; 

Peopyright@ 1984 by Iraklis Kollias": 

eae  *¢ OF ALL C.P. ",N 

: dim BS[N,16],Q$(8],X(N],Y(N],Z(N],A(N,8],S[4],BIN,3],CIN,3],A$ [10,60] 
megsb "INSERT" 

Seps1 DIGI." 

: gsb "RECORD" 

dsp "FIRST STEP DONE ++ USE FILE #l1";end 


’ PILE QR KKKK KK KK KK KKK KKK NN 


WowtDu & WwW POF © 
ee ee ee e@ 6 eo se 


mos "INSERT" :f£xd 0 

ms & “insert name of C.P." 

Mmereror I=1. to Nrprt;lsent, "NAME OF C.P.+ (max # of char=16)> + + +",B$[- 
eas next I 

14; ¢ "Insert the measured coordinates" 

mor iL=1,.to N-fmt 5,f£2.0,".",cz;wrt 16.5,1,85[1]} 

moe ent “Xs>"*,X[I];ent "Yeo" ,Y(I]sent "Z+>", Z[1) 

Mement ,.£2.0.,x,cl6,3£9.2,z:;fmt 2,czwrt 7.1,1,BS{[1),X{1) ,¥{1),2(1) 
meeent “ERROR ? (Yes or No)",Q$[1] 

Rei can (OS fl,1))="Y"swrt .7.2, ERR" s;wtb 7%, 10,133gto -3 

eet 7 ,10,13:next I 

me: % “Output and corrections if any" 

meeiwtb 7,12:fmt 1,20x,52"=" 

mee fmt 2,20x,” CONTROL POINTS MEASURED (X,Y,2) COORDs" 

24: fmt 3,45x,"Xcoord",4x,"Ycoord",4x,"Zcoord" 

Meecme 4 ,21x,£250,".",x,cl6 3£10.2 

mos wrt 7elewtb 7,1l0;wrt 7.2;wtb 7,10¢wrt lcedse Wee LU 

merwrot 7.3e-wtb 7,10 

28: for I=] to Njwrt 7.4,1,BS(1],X[1],Y(1],2{I}¢next 1 

moe wtb 7,l0swrt 7.liwtb 7,12 

30: ent "ANY CHANGE ? (Yes or No)",Q$;if cap(Q$[1,1])="Y";sfg ligsb "ERR 


Muetoufigils;cfig l;gto -5 


B2:; ret 
33; 
34; "BERR": 


35: enp "ROW #",rO0;enp "COLUMN + Ard 
mo: if rl=lsenp X(r0] 
mee if rl=2;enp Y[r0} 
me: if rl=3;enp Z[r0]) 
89: ent "MORE CHANGES 
40: ret 

4l; 
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eee DIGr”:prt “USE THE ““tree* ORIGIN" :spc 2 

mes fOr T=l1 to 2 

eeeidt I=l-eprt “FIRST FILM..." 

Set PL =2sDrt “SECOND FILM..." 

mos, ds “digi TOP END OF FILM";red 4,rl,r2;beep 

@7: dsp “digi BOTTOM END OF FILM":red 4,r3,r4;beep 

mee dsp “digi LEFT LEFT END OF FILM";red 4,r5,r6;beep 

me: asp “digi RIGHT END OF FILM";red 4,r7,r8;:beep 

Mery f T=) er2-(r2-r4) /2+S([(1} sr 7-(r7-1r5)/2-5 [2] 

Mmeeert [=2 -r2- (xr 2-14) /2>S(3] sx7-(r7-r5) /2>5 [4] 

me prt “DIGIT 3 FRAMES OR #3 ®TIMES THE SAMUEL FRAME “spcieZ 

me tOr J=l to 3 

Ba; for K=]l to N 

medcp ™ DIGITIZE ~",BS[(K];red 4,B[K,J},C{K/J] ;beep 

396: next K 

emer "PR R~O'R ? (Yes or No)", OS; if ¢cap(O$S[1,1] )=H="Y"szgto —-3 
Bos next J 

meee "A{L,1) =Ucameral, A{L,3})=Vcameral” 

meee “A(L,5] =Ucamera2, A[L,7] =Vcam2era2" 

Meee ATL, 21 ,A(574} 7X(b,6) ,AlLL:8] = S.D.” 

mee fOr L=l to N 

Merch) MEAN’ (B(L;l1);B0L, 2) BIL, 3],A(L,41-3] ) 

meet “MEAN (C[L;1],C{b,2],CLU, 3l,Altb,41-FF) 

Mec) “Sp (BIL, 1] ,B8(b,2],BIL,;3),A(L,41-3),AI[L,41-2])} 

Mec) “Sp 2(CclL,1] ,C(L,2},Cfh,3)] AlL,41-1] ,Al(L,41}) 

oy: next L eee f : 
68: fmt 2,c9," CONTROL POINTS DIGITIZED COORDs (in, digi-uniits,) 
69: fmt 5,x,f2.0,".",Cl17,£9 .2,£9.4,£9.2,£9.4 

mi: gsb “PRINT" 

71: for Y=l to 100;dsp "second filn now";beepswait l0;next Y;next I 


Be: ret 


74: "MEAN" : (pl+p2+p3) /3+p4;ret t 
"Sb": {(( (p4—-pl)* 2+ (p4-p2)“2+(p4-p3) 2) /2)>porret 
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m7: "PRINT": 

fee 1, 58"="7tmt 3,.24x,."Xcoord”,4x,"S.D.",5x,"Ycoord”,4x,."S.D." 
Re OX 1g pf 20; a 7p aXe Op phe, aX, lode Eo, lax, Oe ten 
meeeeweb 7,27,69;for A=1 to 20;wtb 7,32:next A:wtb 7,27,77 

Seeewrt 7.lswtb 7,10 

Meer L=lewrt 7.2," FIRST-" 

Mee dt L=2swrt 7.2, "SECOND-" 

meee 7,lOswrt 7.lswtb 7,10swrt 7.3-ewrt 7.4,41-3,41-2,41-1,41;wtb 7,10 
mee fOr J=l to N 

ee fe, J, boi J)|,A{J,41-3] ,AldJ,41-2].,,AlJ,41-1j} ,A[J;41) 

me. next J 

meewebD 7,l0swrt 7.lswtb 7,12;ret 


90: "RECORD": 14N+62r0;24N>rl ;64N>r2 

Meet STORAGE SPACE” «prt, "NAME OF (C.P.",rO0;prt "X,Y, 2 coord”, rl 
eo bCis.. COOrd’, (2; prt. ----4-- + ‘ 

Meort "Mark files if you haven’ t done it yet!!!" 

Mme. prt "---------------- i 

mmeort “T¢1—-F#O (1.D.info.)",659 

mos prt “T#1-F#1l (C.P.names)",r0 

Ree rt “TY1—-F#2 (X,Y,zZ coords) ",rl 

Mee prt "“TH1-F#3 (U,Vcoords)",r2 

mes prt “THI-F#4 Us,Vscoords",208 

100: for R=l1 to 10sent "I.D. INFORMATION(10*60char)",AS[R];next R 
101: ent "ARE FILES MARKED ?? (Yes or No)",Q5 

m2: if ca Slee Ys Ggto: +2 

103: Ree 1,6503;mrk 1,rO;mrk 1,rl;mrk 1,r2;mrk 1,208 
mee: trk l;fdf O;:rcf 0,AS$(*] 

mos trk l:fdf l:rcf 1,BSI*] 

Meeetrk l:fdf 2;rcf 2,X[*],Y([*),21* ] 

Mey: trk l;fdf 3;rcf 3,A{*] 

Meee trk l-fdf 4:rcf 4,S([*] 

109: ret 
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ee ee ee ee ere ee ee ee eee 
CALIBRATION COEFFICIENTS (NO REFINEMENT) +TEST 

Taper DeLatalek. 

track #0 


file #1 


DESCRIPTION 

This program calculates the basic DLT coefficients. 
After the preliminary tests the basic DLT coefficients 
together with the coordinates of the point of symmetry are 
SrOred Win Lite, £4. frack, Ml; -OlmD weal. l.K, tape... Lt Can .auso 
be used for testing the CP's coordinates to determine 
whether or not they are correct. After the calculation of 
the coefficients the X, Y, Z coordinates of the CP are 
Simulated and printed. 
Another test that this program is featured to run is the 
Simulation.of the X, Y, Z coordinates of the two cameras 


with respect to the tree origin. 


Several runs with all the points and with less points are 
advisable to ensure that no errored CP participate in the 
Galetiation of the coefficients. 
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° % x 
PCAL 
e OD 
; dim 

ent 

dim 

dim 

2 ent 

ent 

ed im 

mes trk 
mo: trk 
mo: trk 
mo: if 

ma: if 

mo: ara 
mo: "CO 
mr: for 
me: 2(L 
mo: for 
: WxX[ 
21: 0+B 
me: —-A[ 
eo: -A[ 
24: 1+B 
o>: nex 
mo: trn 
my: for 
ao: G[l 
Bo: G([5 
30: G[9 
me: rl/ 
B2: nex 
mo. fmt 
Ba: fmt 
B>: fmt 
36: fmt 
mo: fmt 
Bo: fmt 
B92: wrt 
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eget” ASD Lhe. k TRK Oy FILE LR AEKRKKKEKKK KK KK KK KKK KKM 
IBRATION COEFFICIENTS (NO REFINEMENT) +TESr": 
yright@ 1984 by Iraklis Kollias"; 
A$ [50]:" CALIBRATION WITH DLT (NO REFINEMENT) “AS 
"% OF ALL C.P.",r99:r99A 
Mie YtAl, 2{A],A(A,8),H(4,3],113,41,N[3,3],314],K13),0(21 
BS [A,16],0${16]} 
Piogtl ilk VERY BIRST RUN? .(Y or N)",OS 
"¢ OF CONTROL POINTS TO BE USEDes 
Pees hl AoC le chZGheG bLl), (11,12), F{ 11), uo, lis bene eal 
Paar ol, BS 
Bey Lateee 1X | Xho 1% Je 2 L* | 
Meee Lis oS gay [ *] 
PaO oi, bia) = Y tty ie lLdt 4, S.[%) 
Seca lin) te tlk leddat 4,L(*},S1*] 
S*P;.01-W 
EFF" :0+r0 
=a tO, dy DV, a2 
-1)+1+A;A+2+B;1+tr0+r0 
Pe ie On Cae ) tC J. 
Pier Cit pl] +B (9,, 5.) sWY LI) +Bi 1, 2)>Bild 761 eWw2 oe oe eto 
[1,5]+B[(1,6]>B[I,7]+B(1,8]+B[J,1]*B(J,2)+B[J,3)]+3[J,4] 
eK) +B, oe A lol AWAY dl eB, LO) s—AlT AIWZ ft) sei) 
I, B) WX[(1I]*B[J,9] ;-A[I, BJ) WY[I]+B[J,10];-A[1,B]WZ[I] +B[J,11] 
[1 4)+815,8.) ;Ab1,Al>C(E):AlL, 8) +Ci Jj) 
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ate Ll sis LO tahdaled Oya emex't oJ 
1G FGI21G (LO bGI3)G[ Lilart 
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me -C.P.RECOVERING": 

mee fOr I=] to 20;0-rI;wtb 7,32:next Me WED 2d, 2) pad 
meerimt 1,£2.0,z; ake epLoepa. fat or 74"- ss ene Cd Bae Ss 
Seeeeiat 5 ,18x,c:; ene GO; 968 

meeewret 7.4;wrt 7.5, AS; WECaL 3 

49: wrt 7.6, "Xm", "Xs", "dX", "Ym", "Ys", "4¥", "zm", "25", "42" ewrt Tha 
me: WtD 7,27,86,int (12/64) ,int (12) 

meer or I=1 to r99:wrt 7.1,I1 

mee LOor-T=1- to 3 

mee b(1,T+8j)A{1,1]-L{1,T)+H[1,T] 

Be: Lil, T+38JAl[1,3]-L[{1,T+4]+H[(2,7T] 

Meee ic, (+SJAlI,5)-L[2,T] -H{(3,T) 

ees L(2,T+S8jJAl[I,7)-L[2,T+4]-H [4,7] 

(Mees next T 

Meret 1,4)-A[1,1}+J7[1]);L[1,8}]-A[I,3])+J[2] 

Wee (2,4) -A[1,5)*J[(3]);L[2,8]-A[I,7] +3[4] 

60: trn HeI;mat Il+N;mat IJ+Csinv N>N;mat NO+K 

61: XA{I}+rl;K{[{1] /Wer2;rl-r2-r3 

62: Y[I)+r4;K[2] /Wer5;r4-r5>r6 

63: Z[1])+r7;K[3] /wer8sr7-r8>r9 

64: if I<=C;r372+r10erl0;r6° 2+rlierll;r9° 2+rl2+r1l2 
Mos if I>C;r3°2+r1l3erl3;r6°2+rl4+rl4;r9° 2+rl5-r15 
Meeeetor Y=! to 9swrt 7.2,rY;next Y 

me wtb 7,10,13;next 1 

Meewr eu, .Setmt bebo pfllhe2p2£24 62 

eos rl0O+rl3*rl6;:rllt+rl4*rl7;rl2+ri5+r18 

mre ¥(ebO/C)>rl0:7(r11/C)srll;7 (r12/C)2rl2 

mee, if (C<r 997 (r13/ (r99=C) )er1337¥ (c14/(699-C) )>r14; ¥ (r15/ (r99-C) )>r15 
2: ¥(rl6/r99)erl6;7V (r17/r99)erl7;¥ (r18/r99)+r18 

meee wthbo?, 27786, int (16/64) ,int(16) 


mee SO<r99e8 wrt’ 7.1,"RMS OF C.P. ",xl0,rli,ri2 
Moe if C<r99ewrt 7.1,"RMS OF UNKNOWNS",1r13,r14,r15 
me wrtosvsl,"TOTAL RMS MF GAG pap 


97: fmt 1l,c:wrt 7.1,"*** (Measurements in cm.) ";wrt 7.4 
me: webo7Si2Fwtb 7,27,69 

ms ent "SHALL COEBFF.. BE RECORDED? (Y,N)",O9_ 
m0: if cap(9S[(1,1])#"Y¥";dsp "OH!! THANK S"“;wait 2000;gto +4 
Pe orte*Usl*vsi;0s2,Vs2" 

m2: for Y=l to 4:prt P{Y],S(¥]lzspc +PIY]+Si¥];next ¥ 
Mercer K 1 cert 4¢Li*),sd*4 

B4: trk 0 

85: %$ "CAMERAS COORDINATES" 

Meet or  [=1,toi2-for J=lvto 3 fee 

87: L[I1,J]*N[1,J3] sL[1,d+8] *N(2,5] ;L{1,5+4] *N[3,J] 

BS: next J 

9: -L{[1,4]+O[1] ;-120[2);-L[1,8] 7913] 

90: inv N*N;mat NOK 

al: Pete COORDS. Of Sree Hi" 71 

92: prt K({l]/w,K[2]/wW, K[3]/W 
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SSS 
REFINEMENT WITH Karara's MODELS (main program) 
Tape: Deb. Th oK. 

track #0 


Evie “#2 


DESCRIPTION 

This program was designed for the calculation of the 
DLT coefficients whenever one of the mathematical models for 
image refinement is used. 
User enters 
The # of all the CP in the data file 
The Model # (2,3,4,5,6) to be used 
The # of CP to be used for the calculation of the 
coefficients 
The computer 
loads the CP's data 
loads the appropriate model to be used 
calculates the coefficients 
simulates the X,Y,Z coordrnates of all the CP 
prints “in the small printer the sum of the residuals after 
each iteration 
prints the measured, simulated and their differences 
prints the RMS error of the CP 
prints the RMS error of the unknowns 1f any 
records the calibration coefficients in track #1, file 
Balm). 86 (rir) * 47 GLv).,' BV), oC) 


After the tests of the points are finished it is adviseble 
to run the program with all the correct points in order to 
calculate the coefficients that will be used in the recovery 


of the coordinates of points of interest. 
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"REFINEMENT WITi Karara’s MODELS": force 
Mm COpyright@ 1984 by Iraklis Kollias": 
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SUBROUTINE MODEL II 


Tape: Deb... 1y,.Ke 
track #0 


File #3 


DESCRIPTION 
This subroutine is loaded at the end of the program in 


file #2. The basic DLT equations are used together with 
equations (A.49) and (A.50) for the calculation of the 
calibration coefficients. 
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SUBROUTINE MODEL III 

Tape rDe le reek: 

track #0 


file #4 


DESCRIPTION 
This subroutine is loaded at the end of the program in 


file #2. The basic DLT equations are used together with 
equations (A.51) and (A.52) for the calculation of the 


calibration coefficients. 
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SUBROUTINE MODEL IV 

Tape Diauerohek: 

track #0 


file #5 


DESCRIPTION 
This subroutine is loaded at the end of the program in 


file #2. The basic DLT equations are used together with 
equations (A.53) and (A.54) for the calculation of the 


calibration coefficients. 
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SUBROUTINE MODEL V 

Tape, DAL. Td .K. 

track #0 


file #6 


DESCRIPTION . 
This subroutine is loaded at the end of the program in 


file #2. The basic DLT equations are used together with 
equations (A.55) and (A.56) for the calculation of the 
calibration coefficients. 
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SUBROUTINE MODEL VI 
Tape Deli.t. lak. 
track #0 


Eileuf7 


DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #2. The basic DLT equations are used together with 
equations. {(A.57) and (A.58) for the calculation of the 
calibration coefficients. 
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MODEL": “REPINEMENT WITH MODEL #6 (KARARA) "+aAs 
“Copyright@ 1984 by Iraklis xKollias": 

ror i=l to 2;0*r24>r 43+r25+r26:10* 6>r90 

pei" :beep;fxd O;l+r24+r24:dsp " PLEASE DON’T TOUCH>",r24 
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Tape Di L.Tal ok. 
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DESCRIPTION 

With this program, the user digitizes the points of 
interest from films from both cameras. 
The maximum # of frames for 20 points in each frame is N = 
87 frames when the memory capacity of the computer is 64 
Kbytes. If more frames are to be analysed the digitizing 
Should be done in groups of points. The latest was the case 
for the hammer throwing analysis where 222 farmes were 
analyzed for subject 2. 
User enters 
# of analyzed points from each frame 
# of frames 
5 lines of comments with 60 characters each 
the names of the points for analysis 
track and file to record the U,V coordinates of the points 
from) thes filmoof*thedfirst camera 
track and file to record the U,V coordinates of the points 
from the film of the second camera 
marks the files for the coordinates (the # of bites is 
printed) 
The computer prints the comments, the names of the points 
and the files where the data are recorded. 
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The film from the first camera is to be digitized first. The 


film from the second camera is to be digitized second. 
eOUIE Cre be 0 A as ee eee 


a a i a ne ee ; al 
hss f 7 all j aah, 


Pa ' { rm 
a 


Tad A0F 1010 ATAT ee 
M.1/T.d,0 ager | Abs Dibra 
0% Hosts a | +imal “a ee 
8% efi? 


ee ey te 


10 einiedg ede eoulsipibs 19au arr ams Py 
aaemed> tod ov ; iene 
- Mai emes® Hoge ni esa OS 307 to, \m ms 3 


63 al teduqmpol eft to ytiosqes yz It oedy gem 
onizisiphs atid bhevians ais a4 e768 mio ahs 
s2s2 eff eaw Seedal sgT .asniog Yo aque i’ = i - 

2184 eanr5? S85 sede ae snd eee is 


yy’ , te 


rig 


rt. ee 
es ' 


od? .te3ii bssizipib sd o2 
Fruseen besitipth ed 


en 


: 415% ' ~~. i ; Ci. te 
’ 7 “We AY 
te iV Wioae 


PwWoWYHS UYU &wWwHPY © 
ee ee ce@ eve 


Cue? 66 ce. ee 68 «Use 


B3: 


22408 


% NkkKkK* REKKEKKEKKEKEKD TT. bigenrees , TRK 


me OATA DIGI FOR DLT" 

meooyright@ 1984 by Iraklis Kollias”: 
ent "# OF ANALYSED POINTS IN each ERS 2G 
ent ~% OF- FRAMES. [20points. max#=87]",N 
fe 0sdim A{[N,G,2],D{N,G,2],AS[G, 30],Q0$8{20],8$[(5,60] 
Bom I=1 to 5sent "1I.D. jand COMMENTS", BS [tj] mext J 
meee, 2 UX, CsfOr I=] to Siwrt 7.2,BS[1I] next I 

mem tO Grport Is:ent “NAME OF POINT" ,AS|[1I)*nexteal 
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Per f=1 to 2 

for J=l1 to N 

mOr K=1l to G 

dsp “DIGITIZE POINT",AS [K] ;red 4,U,V;beep 
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mee t=2°UsD(3,K,1] :V>+D[J,K,2] 

mext K 

Otel ERROR ?? Yes: Or ...)", O05 

Pe tcaplos t,t) )="Y";dsp "DO IT AGAIN" ;wait I500;gto -6 
Peo. NE UMRAME NES) 0" ,Jitlizwait, 1500 snext: J 

mom lmeto) 0 eds pi “LOAD «THE SECOND F LUM” pwa dt 555¢next, K 
next I 

GsounNINSERT UV TAPB; © CONTINUE ";sto 

prt "EACH FILE HAS",N*G*16,"BITES";spe 2 

Gap femekkK THE FILES then “CONTINUE “sstp 

en “trktBORMU,V from lst CAMERA", r0;ent “FILE, rl 
Ee tas ar Stier 1. Al [i* | 

ent trk for U,V form 2st ‘CAMBRA",r2;ent "FILE",r3 
mee tra) ty Gticr 3,/ [* } 


’ FILE S8®®®XR KK KKK KKRKEKKN 


Pneee20x , "U,V FROM FIRST CAMERA IN PoRACK “toy fo 2 Up ete uae pote 


mee kat, wOYr L 


fmt 2,20x,"U,V FROM SECOND CAMERA IN TRACK#",£3.0, "FILE #",£3.0 


inl) Acie 2 03 
wtp) tll dsp e’ DONE" end 
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X,Y,2Z,RECOVERY WITH DLT (KARARA's MODELS ) 
Tape ;-UieLinslain Kes 

track #0 

File #9 
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DESCRIPTION 

This program is for the recovery of the coordinates of 
the points for analysis. The basic DLT equations with the 
appropriate model for image refinement are used for this 
purpose. User. enters 
model to be used (2,3,4,5,6) # of analyzed points in each 
frame 
# of frames 
track and file where the U,V coordinates from the first 


camera ' 
track and file where the U,V coordinates from the second 


camera 
track and file where the X,Y,Z coordinates of each point to 
Demotored.. tail paint Wlid. pe. Pecorved snore 7 ise, that 1s 
Miebie ialrol tile Wild De, POInk mi Cleemignt eves inate 
second file will be point #2 (ie left ear)..... in the last 
file will be point #last (i.e. coca cola label of the left 
shoe ) 

The computer 

loads the coefficients of the DLT é 

loads the subroutine with the appropriate model 

loads the U,V from both cameras 

calculates the X,Y,Z coordinates of each point 

marks the files for storing the X,Y,2 coordinates 
rPecordssthesx,Y,2 Coordinates 2n files 


The X,Y,Z coordinates of each point are recorded in its file 
with an array of Gl3, N] saad N = # of frames and G[1, N] = 
ReaGh2) Ni] «oY and G[3, N] = 
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% NEKEKKKEKKEKKKKKKK KK | Te Tipe Pat Kise; TRK Ore FILE QGRKEKKKKEKKAKKKKK KK KN 


feA, X76, RECOVERY WITH’ DLT (KARARA’s MODE [soins 
feopyrighta VSc4 by lfaklis Koll ias™. 

ment MODEL TO BE USED",M;12+C 

: $ “Loads the epee tvcienta! 

if M=3;14+C 

>: if M=4;16C 

if M=5;18-+C 

if M=6;22+C 

emer ie, Cc) ,L, 2,11} ,S [4] 
Meeecrk il;fdf 4;:ldf 4,L[*],S 
Merecrk l:fdf- 3+M;ldf 3+M,P[ 
memoin 1(4,3},1[(3,4],N(3,3] 
me: ¢ "Loads the MODEL" 

mee trk O;fdf M+7;1ldf M+7,39 
eee PP sdsp, "INSERT DATA TAPE; “CONTINUE "stp 
16: ¢ "Loads the data" 

me ent "# OF ANALYSED POINTS IN each FR." 
Ic: ent "# OF FRAMES [20points max#=37F]" 
Be: Ge+rO0;sdim G[3,N],AI[N,G,2],D(N,G,2] 
_mormce tr kK with U,V from ist CAMERA",rOQsent “FILE™,cletrk rosiar rial) 
Mmment “trk with U,V form 2st CAMERA",r0zent “FILE” ,rl;trk ru;lat r2, D|*| 
Meeeent "TRK# TO RECORD Xi,Yi,2i",r51 

Ment "PIbs ¢° TO BEGIN WITH", r52;r52-l+r52 . ; 

Memmeore "STARTING FILE=",r52+l;spce ;prt “# OF FILE=",6 

me prt “BEACH FILS HAS" ,N* 24, "BITES" ;spc 2 

Mme dsp “INSERT X,/Y,Z, TAPE; “CONTINUE © pers Ge) 

mm: for A=1 to r0;for I=1 to Njdsp A,I 

ms A{I,A,1})>M[{1,1] ;A[1,A,2]*M11,2] 

meee Di 1,A,1)+M[(2,1) -D[1,A,2}>M[2,2] 

ee M{1,1)=S[1])>M[(1,3] 2M(1,2]-S(2)-4[1,4] 

Me 1(2,1}-S[3)}om(2,3) pM(2,2)-S(4]*M(2,4) | 

92: for K=1 to 2;7 (M[K,3])~2+M[K,4] 2)*M[K,5] ;next K 

* ash ’RCv" 

4: next I 

SB: if At+r52=0:rew:trk r5l;mrk 1,N*24;rcf Atr52,G6(*] 

Bo: if At+r52#0;trk r5l1;fdf Atr52;mrk iE N*24srcf Atr52,6G[* } 

m7: next A 

38: asp "DONE"send 
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SUBROUTINE MODEL II 

Tape D.&.T. EK. 

track #0 


file #10 


—_—_ nee SO ee ee 


DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #9. The solution is based on equation (A.79). The basic 
DLT equations are used together with equations (A.49) and 
(A.50) for the calculation of the 3-p «coordinates. 


v) % “MODEL #2" 

l: gto tpi 

2°: PRC Vs 

Beeror J=)l to 2:J3+2+L 

mercor K=l) to G:P[y,K)>rkKsnext K : 
5: A[J,1]+O;M[J,2)+P:M[J,3]7+O:;M[J,4] >ReM[J,5] 5 
6: O+rl2RS°2>r20; Ptr1l2RS° 2er21 

97> rl—-r9r20-[J,1) ;r5-r9r21-HI[L,1] 

6s r2-rlOr20-0[(J,2] ;r6-rl0r21>+H(L,2] 

9: r3—-rlir20eH (J,3]):r7-rlir21-HI[L,3] 

iO r20-r4>J [J] 21-r8+7 [Lb] 

ms next “J ‘ 

12: trn Hel;mat Ill-N;mat IJ*O;inv N+N;mat NO+K 
13; for T=] to 3;K[T)+G(T,I)snext T 
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SUBROUTINE MODEL III 
Tape? D),l:.T.,l.Re 
track #0 

file #11 


DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #9. The solution is based on equation (A.79). The basic 
DLT equations are used together with equations (A.51) and 
(A.52) for the calculation of the 3-p coordinates. 


MODEL #3" 

Wy) S7 0) 
Viens 

J=1 to 2;J+2-L 

K=]1 to C;P({J, K] erkK;next K 

pL )sO7M(3,2)>P2M[(JI,3170;M[J, 4] *RIM[J/SJ+S 
(126 724r13S44+r14S8" 6) >r20; P+ Riel 2S i2+r3se44rl 4s 6)>r 2) 
r9r20+H[J,1) ;r5-r9r21l+H(L,1] 
rl0r202H (J,2] ;r6-rl0r21+H[L,2] 
rllr20+H[J,3];r7-rllr21-HI[L,3] 
O-r4+J [J] ;r21-r8-+J [L] 
xT iJ 
n H*I;mat IH*+N;mat IJ+O;inv N+N;mat NO+K 
@ T=) to 3-K{(T)*GiT, 1d) ;next 7 

c 
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SUBROUTINE MODEL IV 

Tape D. i, Ti kK: 

track #0 

file #12 
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DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #9. The solution is based on equation (A.79). The basic 
DLT equations are used together with equations (A.53) and 
VAtod) thor ithe calculation of the 3-p coordinates. : 


ms "MODEL $4" 

l gto "pp" 

me "RCV": 

me for J=l to 2:J+2-L 

®: for K=l1 to C;P[J,K]+rK;next K 

5: M[J,1]+O;M[J,2]+P;4[J,3]>Q;4[J,4]*RiM[J,5]45 

6: O+0(r12572+r133~4t+r14876) +415 (S°2+2Q 2) t2r16QRer20 
Me PER(r125 ~2+rl13S 4+r14S 6) +rl6(S° 2+2P 2) t2r150R-r 21 
8: rl-r9r20-N[J,1] ;r5-r9r21+l[(L,1] 

9: r2-rl0r20+H [0,2] ;r6-rlOr21+li(L,2] 

WO: r3-rllr20-H[J,3];r7-rlir21-HI[L, 3] 

Mls r20-r4+J [J] ;r21-r8-+J(L);next J 

12: trn HeI;mat IlfleN;mat IJ+O;inv N*N;mat Nok 

me for T=1 to 3;K(T]+G(T,1I] ;next T 

14: ret 


227), 


rs 


o Oe eee eae neem atime 


VI J2BaoM RUS TORS 
.4,1,.7.2.¢0 eaat 
7% woet? 


STR ef53 


at maypexg 929! 
rived oA. .(3¥s. 


 aeat ) (fG8) 


viv ceaaae Mune ats. 
ome ci Estas 2 Bt 


dn 


al 
i 


‘a) i i : sd al 
ie . _— : bf 


i Ae 
| aT 
= pi, U 
'' : i 


aes ne 


BPRS 


SUBROUTINE MODEL V 
Tape, Dy..b...T. li. Ke. 
track #0 


File #13 


DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #9. The solution is based on equation (A.79). The basic 
DLT equations are used together with equations (A.55) and 
(A.56) for the calculation of the 3-D coordinates. 


O: % "MODEL #5" 

1 gto "pp" 

2: wRCV ts 

Be f£or J=)l to 2: J+2>L 

Meetor K=l-.to C:P[J,K]>rk :next K 

6B: M[J,1J)+O;M[J,2]+P;M[J,3]+Q:M[J,4]*ReMlJ-5]75 . 

6 Oty (r12572+r13S~3+r14S “4+r15S°5+r16S 6)trl7(S 242) 2) +2r18IR+r20 
i: P+R(r12$*2+r13S° 3+r1435 “4t+r15S 5+rl6s 6h trie VS 2tZR 2)et2r CRs re 
8: rl-r9r20+H[J,1] ;r5-r9r21+HI[L,1] 

9: r2-rlOr20eH [J,2];r6-rl0r21l+H[L,2] 

10: r3-rllr20+H[J, 3] ;r7?-rllr21-H[L,3] 

el: r20-r4+J[(J] ;r21-r8+J(L] 

B2: next J 

13: trn HeI;mat IH+N;mat IJ*+O;inv N*N;mat NO>K 

14; for T=l to 3;K(T]*GI(T,I];next T 

5: ret 
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SUBROUTINE MODEL VI 

Tapes D ais. Wale ha 

track #0 


file #14 


DESCRIPTION 

This subroutine is loaded at the end of the program in 
file #9. The solution is based on equation (A.79). The basic 
DLT equations are used together with equations (A.57) and 
(A.58) for the calculation of the 3-D coordinates. 


& MODEL #6" 

gto "upp" 

eRCV": 

Por J=l1 to 23:J+2+L 

mecor K=l to C-P[J,kK]>=rkK next K 
M(J,1)+O;M[J,2]>P:M[J,3}]*Q:M[J,4]*>RsM[J,5]>S 

OF0 (vr 125 -24+r13373+r145 44+r15S°54+r16S 6) 4r17(S 2420 2) 4+2r160Rs141 

Petri 2S 247138 34+r148 “4471 5S" 5+rL6S O)4ri1e(S 24282) 42rl 7 ore? 

> r414+r190° 2+r20P* 2>r39;r42+r210° 2t+r22P 2+r40 

> rl-r9r39eH(J,1] ;r5-r9r 40-1 (L,1] 

> r2-rl0r39-H[J,2];r6-rl0r40+H[L, 2] 

r3-rllr39eH (3,3) ;r7-rl1r40-HI[L, 3] 

r39-r4-J [J] ;r40-r8-J[L] 

next J 

trn H+I:mat IH*N;mat IJ*+O;inv N*N;mat NO+K 

Pormr=1 to 3: K(T1]e+oiT,1)snext it 

ret 
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APPENDIX B 


Moments and Products of Inertia 
Angular momentum 
and 


Instantaneous Angular Velocity 
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This appendix deals with the calculation of the moments 
of inertia of the system athlete plus hammer about a 
Space-fixed orthogonal system of axes with origin at the 


center of mass of the system. 


Moments and Products of Inertia with respect to CMs 

Figure B.1 1S a representation of a jth body segment in 
Space. Let O be the center of mass of the system (CMs) and 
let C be the center of mass of the ith segment. Let X, Y, Z 
be right handed orthogonal system of axes fixed in Space and 
having origin at the point ©. They are translating as the 
center of mass of the system translates. Axes X', Y', Z' are 
Pavallel axes to the K, Y, Z With origin at point C;. 7Axes 
RepyeYcap Zc;;, are the principal axes of the fth segment, 
having origin at the point C. 

The sum of the moments of products of inertia of the 13 
segments about the point O are given by the following 


equations: 


Ix = EEL Cos + ems Conergy ane (Bide 
ly = ELE yn tam pCCgaewh eee ot (axe) 
tee EET 21 + mi (Ci? + Cys?) (ae) 
hey i= Lidds rd tompsaCkneCha 3 (B.4) 
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A ’ 
Pa ZEIT x2 seat pe Ca eaaet (B.5) 


Dye 


3 
= LI" ya tin Cy ean (Bee) 


where: I,, I,, I, are the moments Of inertia of the system 


about the axes XK, Y, Z 

m; is the mass of the ith segment. 

Cyxi, Cyi, Cz; are the coordinates of the center of 
mass of the /th segment relative to the center of 
mass of the system. 


Dis ig) ley A Scyte e ee ee amomerice 
and products of inertia of the jth segment about the 


axes? Kiap Yio nee. 


Moments and Products of Inertia of segments about their CM 


The moments and the products of inertia of the jth 
segment with respect to the axes X', Y', and Z' passing f 


its CM, are calculated using the following equations: 


rom 


Mea lo ey ot 1 yet Oza eee ae (3.7) 
eae, tag OL a2 3s ot Oke CB by 
0 oh i ae: ae og i eatats 2 Fel ee Oa oe (B.9) 
OS = <4? @yo°k xh + @ig 4 “Og lve + O54 *O59°L 24) (B16 } 
Lig = la, eRe 1 © est C.F ass! ee Oa la) Be) 
a —“(aporas acl x Pat 222075 ls yn Miers ee Ogle) aCe te) 


where : I”,,, 1 yir.! zi are the principal) moments. of 
anec hia of the ith segment. These data were taken 
from Whitsett's anthropometric data, (Whitsett, 
1963), and are presented in table 8 in the 


methodology chapter. The principal moments were 


normalized according to the height and mass of each 


athlete. ; 
@11, Mi2r @i3 are the direction cosines of the 


X¥c (ij = 1), Ye (i = 2) and Zc (i = 3) axis relative 
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Bo, the: (Mie SY, #2’ system of axes. 


Calculation of Direction Cosines of the Principal Axes 
The following is a description of the technique used to 
Calculate the direction cosines of the Xc, Ye and Zc axis. 


The vectors representing the axes Xc, Yc, Zc are: 


R20 (ik EO 2 AN CO pak (B03) 
¥c, = A211 + a22) ae Q23k (B.14) 
ZC = @341 Ls Q32]) - a33k (B.15) 


wheres iy oj,nk jane (the junitrvectons ofthe -axes XK", Yip, Z'. 


Direction Cosines of Zc 


The direction cosines of Zc are given by the following 
equations: 


es (B. 16) 
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Direction Cosines of Xc 

Referring to figure B.1, line Xc is the intersection of 
the planes (G) and (E). Plane (E) is normal to the line 
segment Zc at point C. At the jth frame, the equation of the 


plane (E) is: 
Q231X + A@z32y + a332 = 0 CB. 12) 


where: a@31, @32, @33 are the direction cosines of Zc and 

have been calculated previously using Equations 

GRahé)pn(Beltiiand OBete)< 

Plane (G) is defined by the points P(p,,p2,p3), which 

is the proximal end point of the ith segment in the jth 
frame; the point C(c;,c2,c3), which is the center of mass of 
the segment in the jth frame, (c; = cz = c3 = 0); and 
Q(q1,Q42,q3) which was the center of mass of the segment in 
the (jt+/)th frame measured from the point C at the jth 
frame. For the trunk segment point Q is defined by the right 
hip point at the same frame and its coordinates are measured 


from the center of mass of the trunk at the same frame. 


The equation of plane (G) is: 
Ax + By + Cz = 0 (B.20) 


The coefficients of the equations of the plane (G) can be 


found from the following equations: 


A py, +{BipetttG "ps = 0 (Bo 2h) 
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A 94 Cn.8 G25t € qs = 0 CBA 2) 


By setting A = 1 and rewritting equations (B.21) and (B.22) 


in a matrix form, we obtain: 


PZ Ups PB ~PA 


(B23) 
€ 


Gz Q3 mee 
This system can be solved numerically for B and C. Since 
line Xc is the intersection of planes (E) and (G), it can be 
represented analytically by equations (B.19) and (B.20). 
Assume that a point,R(r,;rg,"3) is a point of the line Xc, 
@hen if f, =.q.i, the’coordinates r2.and rx; are given by the 


following equation: 


(B.24) 


Solving (B.24) for rz and r3 the direction cosines of Xc are 


obtained as following: 
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Oi — a eee en (B.26) 
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Direction Cosines of Yc 


The direction cosines of Yc are calculated by utilizing 


the vector product of the vectors Zc and Xc because 
ze 0x Xe = Xe (B.28) 
Using equations (B.13) and (B.15), the vector product of Zc 


and Xc can be written in the following determinant 


formulation: 


2 Sigh & 
Zc x KC = |]@31 32 33 (B.29) 
O11 @12 413 
or 
NCy = (a32°Q13 7 @33°Q12)1 
- (agqea, 30> @33°A11)j 
+) (a31°Q12 -— &32°041)K (B.30) 


From eqation (B.14) and (B.30) we obtain the following: 


ery) ES deere OCs cy Oars Os CBA3) 
eo. =e 3°ARG SU Ase ans (B.32) 
Q5,05 GQEGRQhe =) A525 04% (BY 33) 


The above are the components of the vector representing the 


Yc axis. The direction cosines Yc are calculated by using 


these components. 
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Angular Momentum of the System 


The angular momentum of the System athlete + hammer 


ait 


about its center of mass was calculated with the following 


equations: 


where: Dy, Ly, Ll, are vectors representing the angular 


momentum of the system about its center of mass; 


(B.34) 


(B35) 


(B.36) 


L',j, L'yi, L'x; are vectors representing the angular 


momentum of the ith segment about the center of mass 


of the system, considering the segment as a point 


L’ xi, L yiy Lozi are the vectors representing the 


angular momentum of the jth segment about its center 


of mass. 


Die yay, le (Pare ccalculated sas gto) lows mete rehnn S27 eoe 


the position vectors of the center of mass of the segment 


with respect to the center of mass of the system in two 


Succesive frames (jth and (j+7)th frame) The components of 


the two vectors are the coordinates of the points C and C' 


(figure B.1) measured from the center of mass of the system 


at the (jt+t/)th frame. A vector in the direction of the 
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angular momentum is given by the following equation: 


N = vy; Xx re (B.37) 


Bermnc.., Cy, Creand Cy'y Cy“, Cae be the components of the 


two vectors, then 


No=) ay “ae 4 + aak (B.38) 
Wee >, Gye Cy iCy yp eyC. «Cy 
rely toma OMS Cp eae en Ces Cot 
Gigs @, iC yb ' Cy +Cy! 


The magnitude of this vector is 


N &y(a.2%.+ ape tha, ahGG3 (B.39) 


The direction cosines of the vector are £,, B:, $3 where 6, 


=a, /N,... etc. The magnitude of the angular momentum 


vector 1S: 


fete. mahi cat (B.40) 


where m is the mass of the segment and At the time interval 


between the two frames. 


Finally the angular momentum of the segment about the center 


of mass of the system is: 
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BOES Ae 6 (B.41) 
Lys iseloB. (B.42) 
Lz" = L+8; (B.43) 


The magnitudes of the vectors Eoupabby ,» ubiagareecabeubated 


due to the following equations: 


Ee cS Teiica W x j Seen ag Wyji — ae Gey We j (B.44) 
L yi ‘age Fr Wy j ee Dy ji 7 ie Wo j (B.45) 
L-2i = PESRLOs | 3 ls gO 4 i Bee pee | (B.46) 


where: Tees aon ae Dares and ee a a Dette Us eer are the 
moments and products of inertia of the /th segment 
about the center of mass of the segment; 
O.i, Oyj, M27 are the angular welocities of =the 
segment about a system of axes with origin at the 
center of mass of the segment and parallel to the 
Space fixed axes. 


Be i> Meet CL/ok, Lye ita Bee are calculated as 


it is shownwinvequdtions (B.7) to (B.12). 


Angular Velocity of the Segmental Principal Axes 


Referring to figure B.1, the ith segment moves from the 
wen frame and from its, position Ci xc ave mc) fomchemie tie 
frame and to the position C'(Xc', Yc', Zc'). The direction 
cosines of the axes Xc, Yc, Zc and the axes Xc', Yc', Zc' 
are calculated as in part 3 of this appendix. Let aj;,, aj2, 
a,3; be the direction cosines of the Xc axis, a21, G22, 23 


the direction cosines of the Yc axiS, @3;, @32, @33 the 
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direction cosines of the Zc ass, 96.47 64 oe ba athe 
direction cosine of Xc', B2;, B22, B23 the direction cosines 
of the Yc' axis and B3,, B32, B33 the direction cosines of 
the Zc' axis. The angles 6, of the Xc axis with the Xc' 
axis, 02, of the Yc axis with the Yc' axis and 63 of the Zc 


axis with the Zc' axis can be found by using the dot product 


of the vectors representing the axes and they are as 


follows: 
Gy = Varecos (Gi48 en tsi yet Gy shad (B47) 
Gs =larecos(@. 16 o4 * @250>5 ass Pen) (B.48) 
63 = arccosS(a3;83; + 432832 + @33633) (B.49) 


The magnitudes of the angular velocities of the axes are: 


WM, = On / Re (B.50) 
We = 02/At (Biot) 
Oo Or / At (B52) 


where: At is the time between the Uth and (J+7)th frame. 


The above angular velocities can be represented by 
three vectors each of them normal to the planes formed by 
thevaxes. Xc and Xce\pevo andwYc sacomandeze srespecrively, 
after the they have been translated to a common origin. Let 


Wiiws and ws besthe three vectors, then 
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= X¢ x Xc' (B53) 
RC CG (B.54) 
= "2c x 2c! (B.55) 


(Note: Wig M2, Ws are not necessarily orthogonalmvectors.) 


By using the directions of the axes and the equations 


(B.53) to (B.55) we obtain: 
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(ar2°Bi3 - ae Cae 
CGsG" hia @13°Bii)j 
(dq a" P12 —-@72°Baaqyk (B.56) 


(a22°Ba3 - Q23°*Ba2)i 
(Q21°Bas - Q23°B21)j 
(d21°B22 - @22°B21)k (Be57)) 


(a32°B33 - a33°B32)i 
(a31°B3s3 - @33°B31)j9 
(a44°Bi2 ~ @12°Biidk (B.58) 


Y¥13 represent the components of wi, Yo4, Y22, 


the components of wa and ¥31, Y32, 33 


components of w3. The direction cosines of 


these vectors will be: 
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Mhe total angular velocity vector of the segment about its 


center of mass is: 


Wea Wr is (B.59) 


The components of the velocity vector in the X, Y and Z axis 


are’: 


Wy = @1°O4y7 F Wo°b21 + wW3?S5; (B.60) 
Wy = W1°b42 + W2°b22 + W3°532 (B.61) 
Wr = W1°b13 + W2°h23 + W3°833 (Be 6 2) 


Instantaneous Angular Velocity of the System 

After the moments of inertia, the products of inertia 
and the angular momentum of the system have been calculated, 
the angular velocity of the system about its center of mass, 
can be calculated by solving the following system of 


equations for Q,, Q, and 2,. 


be = ese ae Ty yt Qy Be AE ON (B.63) 
bY Pie fre ty My ey eke (B.64) 
L, = RAIL3G Dy ashy =i Int, oth (B.65) 


Lx, Ly, Lz are the angular momentum of the system; 
I free; 1, and Ixy, Ixz, Lyz ete the moments and 
products of inertia of the system about its center of 


mass; a 
Q,, Qy, N, the requested angular velocities. 
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Figure B.1 Representation of the ith segment in space at the 
jth and (j+/)th frame. 


O(X,Y,Z) = Orthogonal system of axes parallel to the space 
fixed system of axes and having origin at the center of mass 
of the system; 

C(X',y',Z') = Orthogonal system of axes parallel to 0(X,Y,Z) 
with origin at the center of mass of the segment (point C) 


at the jth frame; Rea 
C(Xc,Yc,Zc) = System of the principal axes of the segment at 


the jth frame; ai lys 
CM(Xc' Yo", Zc')e= System of the principal axes ofuthe 


segment at the (jt+/)th frame; 
P and D = The proximal and distal point of the segment. 
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Radius and Center of Curvature 
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This appendix deals with the derivation of the equations for 
the calculation of the radius and center of curvature of the 
Orbit of a point P which moves along a space curve (I), as 
in Figure C.1. The method is based on vector analysis 
(Spiegel, 1959). The kinematics of the point are calculated 


numerically from the time-displacement data of the point. 


Radius of Curvature. 
Let P be a point moving along curve (I). The position 


vector of P with respect to the origin O is: 


me epxieruys! = ozk ee 


Wheres x, y, Zz are the coordinates of 1P; 
i, j, k are the rectangular unit vectors of the 


reference frame. 


Artangent vector. to (I )tat, pompseeas. 


22 = £detey.’)  2ek CCaey) 


Where: x's yl) z’ denote the derivatives of x; y, Zz with 


7 Lf 


respect to time t. 


The arc length s is related to ©F by the following equation: 


dr ds 
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Figure C.1 Representation of the curve [F in the 
space. 


Normal unit vector to the curve; 
Tangent unit vector; 

Binormal unit vector; 

Center of curvature. 
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